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Abstract 

We re-examine the one-loop renormalization group equations (RGEs) for the dimensionful param- 
eters of the minimal supersymmetric Standard Model with broken supersymmetry, allowing for 
arbitrary flavour structure of the soft SUSY breaking (SSB) parameters. We include threshold 
effects by evaluating the /3-functions in a sequence of (non-supersymmetric) effective theories with 
heavy particles decoupled at the scale of their mass. We present the most general form for high 
scale SSB parameters that obtains if we assume that the supersymmetry breaking mechanism does 
not introduce new inter-generational couplings. This form, possibly amended to allow additional 
sources of flavour-violation, serves as a boundary condition for solving the RGEs for the dimension- 
ful MSSM parameters. We then present illustrative examples of numerical solutions to the RGEs. 
We find that in a SUSY GUT with the scale of SUSY scalars split from that of gauginos and higgsi- 
nos, the gaugino mass unification condition may be violated by 0(10%). As another illustration, 
we show that in mSUGRA, the rate for the flavour-violating ti — > cZi decay obtained using the 
complete RGE solution is smaller than that obtained using the commonly-used "single-step" inte- 
gration of the RGEs by a factor 10-25, and so may qualitatively change expectations for topologies 
from top-squark pair production at colliders. Together with the RGEs for dimensionless couplings 
presented in a companion paper, the RGEs in Appendix B of this paper form a complete set of 
one-loop MSSM RGEs that include threshold and flavour-effects necessary for two-loop accuracy. 
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I. INTRODUCTION 



Renormalization group equations (RGEs) have played a vital role for the extraction of 
phenomenological predictions of theories where the physics is specified at a scale much higher 
than the energy scale directly relevant for phenomenology [1]. A very familiar example of 
this is the prediction of the weak mixing angle 9w in grand unified theories (GUTs), widely 
regarded as one of the important successes of the idea of grand unification in the context 
of supersymmetry (SUSY) [2]. RGE methods have become ubiquitous for many studies 
of physics beyond the Standard Model (SM), where relationships between various model 
parameters, renormalized at some very high scale, are specified by the new physics. 

Of interest to us here are supersymmetric models of particle physics [3-7] where super- 
symmetry fixes the dimensionless couplings of superpartners in terms of the corresponding 
gauge and Yukawa couplings of SM particles. Model-dependence arises via the potentially 
very large number of undetermined dimensionful soft SUSY breaking (SSB) parameters that 
reflects our ignorance about how supersymmetry is broken [8]. In practice, one resorts to 
models with underlying assumptions that fix all the SSB masses and couplings, renormalized 
at a high scale, in terms of a handful of parameters [9-12]. These high scale parameters 
have, of course, to be evolved to the low scale relevant for phenomenology in order to sum 
the large logarithms that may otherwise invalidate any (fixed-order) perturbative calculation 
[13]. 

This paper is a follow-up of an earlier paper, hereafter referred to as Paper I [14], where we 
used the seminal papers of Machacek and Vaughn [15] to obtain the one- loop RGEs for the 
dimensionless couplings of the Minimal Supersymmetric Standard Model (MSSM), carefully 
including threshold effects [16, 17] that are essential to include for true two-loop precision. 
Here, we extend the existing literature, and present the one-loop RGEs for the dimensionful 
(in general, complex) parameters including both flavour-mixing effects and threshold correc- 
tions. As discussed in Ref. [14] , these RGEs can then be augmented by the two-loop terms of 
the MSSM RGEs [18-21], without the need to include threshold corrections to these for the 
required precision. Just as in Paper I, where we found that SUSY breaking threshold effects 
cause the gaugino-fermion-sfermion [higgsino-fermion-sfermion] couplings to evolve differ- 
ently from the corresponding gauge [Yukawa] couplings, resulting in an enlarged system of 
RGEs for the dimensionless couplings, we will see that as we decouple heavy super-partners 



2 



and Higgs bosons, a similar situation obtains for the RGEs of dimensionful parameters. For 
instance, the super-potential higgsino mass /i evolves differently from fl, that enters via the 
corresponding Higgs boson mass squared parameter, Also, at scales below the masses 
of the heavier particles in the Higgs boson sector, where just the SM Higgs boson, h, and 
the associated would-be-Goldstone bosons that complete the scalar doublet (see Paper I) 
remain in the theory, we will see that only certain combinations of model parameters can 
be evolved in the corresponding low energy theory. These complications, to our knowledge, 
have not been included in previous studies. 

We use the RGEs for the dimensionful parameters of a general field theory with real spin- 
zero fields and two-component spinor fields given in Ref. [22] to derive our results. In Sec. II 
we first re-cast the results of Ref. [22] in a form suitable for use by phenomenologists more 
used to the formalism with four- component spinors coupled to complex scalar fields and, of 
course, also gauge fields. As explained in Paper I, this re-casting does much of the work 
needed for the derivation of the MSSM RGEs. We discuss particle decoupling in Sec. Ill 
and describe our derivation of the RGEs in Sec. IV. In Sec. V, we numerically analyse the 
solutions of (some of) the RGEs and discuss squark flavour mixing. We apply our results to 
obtain the rate for the flavour-violating decay of the lighter t-squark, ti cZi in Sec. VI. 
We conclude in Sec. VII with general remarks and a summary of our results. The set of 
RGEs that we obtain are listed in Appendix B. 



II. FORMALISM 



Luo, Wang and Xiao [22], that we use as the starting point, write the Lagrangian density 
for a general field theory as. 



£(2) = ti^yo^^i^, 
1 - 

~ 2 . 



1 



1 



P -r -D^(f)aD''(f)a - - 



2! 



(1) 



1. 



- ( ^YlAplCi)q(t)a + h.C. ) - -Thafec0a0fe 



2 - pq^P ^ T ^Tu, . . 21 

where the matrix ( = ia2 is introduced to make the spinor bilinear Lorentz invariant. The 
gauge interactions for the real fields (pa and two-component spinors ipp are contained in 
the corresponding covariant derivatives, the matrices mj and Y"- are symmetric (but not 
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Hermitian) , while the entries of the scalar mass matrix as well as of the trilinear and 
quartic scalar couplings h and A are symmetric in all their indices as well as real. 

The Lagrangian density can also be written in terms of complex scalar fields and 
four-component Dirac and Majorana spinor fields as [14], 



1 

2 L 



Mk 



+ 



riSa6$a$b + h.C. 



+ i {Xl)^^^M,PL^Mk^a + I (KD^^^MjPL^Mk^i 



h.C. 



+ 



2! 



$];H„bc$6'^'c + h.C. 



1 1 
2! 2! 



3! 



Kbcd^i'^b^c^d + h.C. 



(2) 



The Yukawa coupling matrices U^'^ that couple Dirac spinor fields to complex scalars, and 
the matrices Va and that couple these scalar fields to one Majorana and one Dirac 
field, have no particular symmetry (or Hermiticity) properties under interchange of the 
fermion field indices j and k. These indices label the fermion field type (quark, lepton, 
gaugino, higgsino) and also carry information of flavour and other quantum numbers {e.g. 
weak isospin and colour). On the other hand, the matrices X^'^ that couple scalars to two 
Majorana fields are symmetric under j ^ k because of the symmetry properties [4] of the 
Majorana spinor bilinears that appear in (2). The scalar mass squared matrix — notice 
that we have used a different font in (2) to differentiate this mass term for complex scalar 
fields from the corresponding one for real scalars coupled to two-component spinors in (1) 
— is Hermitian in the scalar field indices, a, b, while the matrix B is symmetric. The CP 
even and odd Majorana fermion mass matrices [4] mx and are both Hermitian and 
symmetric. The trilinear and quartic scalar couplings A, A' and H are symmetric under 
interchanges a ^ b and/or c ^ for A, under interchanges of b, c, d for A', and finally under 
6 ^ c for H. As noted in Paper I, while (2) is not the most general form for the Lagrangian 
density, it suffices for the derivation of the RGEs of the MSSM with i?-parity conservation. 
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The evolution of the fermion mass parameters in (1) is determined by the /^-function [22], 



(4vr)' /3„ 



'/ ll-loop 



=^ [Y^(F)m/ + m/Y2(F)] + 2YV}y' + ^Y^Tr 
-3/{C2(F),m^} . 



Y^Vf + m"^,Y' 



(3) 



Here, C2{F) = t^t^ is the quadratic Casimir operator for the fermions, and Y2{F) = Y'''''Y''. 
Notice that we have altered Y|(F) (recall Y2(-F) is a Hermitian matrix) that appears in 
the first term of the corresponding equation in Ref. [22] to Yj{F). This ensures that the (3- 
function for uif has the same symmetry property as mj in (1). We can write the "Yukawa 
couplings" on the right hand side of (3) in terms of the corresponding couplings in the 
four-component notation by introducing 



and replacing Y" with a (3 x 3) block matrix containing U^, U^, V^, W^, and as 
worked out in Paper I. 

The fermion mass matrix can similarly be written (in this same basis) as. 



m/ 







where in the context of the i?-parity-conserving MSSM, all the entries except {mx — ii^'x) 
vanish because gauge invariance precludes the corresponding fermion bilinears. From now 
on, we retain just mx and in our analysis. We can now substitute m/ as well as the 
form of the "Yukawa" matrices from Paper I into (3) to obtain the RGEs for the matrices 



m 



T 
W 



VI/ y 







[my - tniy} 



\ 



J 



(4) 
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mx and m^. We find, 



(4vr) 



dmx 

dt 



and 



(4vr) 



X 



dt 



+ (mx + zm'x) (WfeWj; + Vf V; 
+ (mx - zm'x) (W,*W^ + V^V 



Xfo + x^ 



W^W^ + V^V, + Xl^Xl + X^^X^ ) (mx + tnix] 



[X^ (mx + im'x) X.I + X^ (mx + im'x) X^ 
+X^t (mx - ^m' ) X^t + X^^t (^^ _ ) x^t 



+ 



X.^Tr jx^t (mx - ^m'^) + (mx + ^m'^) X^} 
+ X^Tr {Xf (mx - im'^) + (mx + zm'^) X^} 
+ X^^Tr I (mx + im'^) X^ + Xf (mx - im'^) } 
+Xf Tr I (mx + im'^) X^ + X^^ (mx - im'^ 
- 6/Cf (F)mx , 



(W,W; + V^V,* + Xl^l^ + X^X^t) (mx - tin': 
- (mx + ^m:^) (w,Wt + V^V* + X^X^ 
+ (mx - ^m;^) (W*W,^ + VtV, + X^^X^ + X^^x^) 



X) 



WlWl + Vtv, + X^tx^ + X^tx2 j (mx + ^m^ 
+ 2 [X^ (mx + tm'x) Xl + X^ (mx + ^m'^^) X^ 

-X^t (mx - tm'x) X^^ - Xl^ (mx - ^m' ) X^^' 



Xj 



(5) 



(6) 



+ [xITt \xl^ (mx - ^m:^) + (mx + ^m'^) X^} 
+ X^Tr {X^t (mx - im'^) + (mx + im'^) X^} 
- X.^^Tr |(mx + im'^) Xl + Xf (mx - im'^)} 
-X^tTr |(mx + ^m^) X^ + X^t _ ,^/^) 

where t = \nQ, 02"^ (F) is the quadratic Casimir for the Majorana fermions as defined in 
Ref. [14], and there is a sum over all gauge group factors in the final term of both (5) and 
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(6). Note also that in the MSSM with i?-parity conservation, the trace terms in (5) and (6) 
vanish. This is because, as seen in Eq. (2), X only connects higgsinos to gauginos, while 
m^^ never connects higgsinos to gauginos. Therefore, a trace of the product of these two 
matrices is always zero. 

Turning to the RGE for the trilinear scalar couplings, we first write the one-loop RGE 
for these couplings in (1) as [22], 

(4vr)' /3h.,Ji„,„,, = AL - ^nabc + AL - ^g^Kbc , (7) 

with 

Aafec = 2 ^ Aafee/he/c , (8) 

perms 



Uabc =\ Tr{m;YtY,Yt + Y.m} Y,Yt} , (9) 

perms 

AL =\ [Tr {Yt,Y, + YlY,,} h,,^^ + Tr {Yt,Y, + Y+Y,,} h,,, 
+Tr{Y^,Y, + YtY,,}h„,,, 



(10) 



Kbc = [C2{a) + C,{h) + C2{c)] Kbc , (11) 

where a', 6', c', and the indices e and /, are summed over all the scalars in the theory. 
Also, there is an implied sum over all gauge group factors in the final term of (7), and 
C2{S) = t^t^ is the quadratic Casimir operator for the scalars. All this is exactly as in 
Ref. [22] except that we have written the expression in (10) in a more transparent form. We 
note that the /5-function in (7) is manifestly symmetric under any permutation of the three 
indices. 

In order to convert this into our complex, four-component notation, we first expand out 
the trilinear scalar term in (2) in terms of real spin-zero fields {i.e. the real and imaginary 
parts of our complex fields), and compare this with (1) in order to obtain iiabc in terms of 

Habc = -^-^ [~^aBbRCR + ^aRbjcj + ^^aRbRCj + ^^aRbjCR] ■ (12) 

We can now write down the RGE for iiabc using the same relations for the two-component 
Yukawa matrices as before, and substituting the trilinear and quartic scalar couplings habc 
and Xabcd in terms of the corresponding couplings Habc and Aabcd or A^^^^ on the right hand 
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side of (7), taking care to sum over real and imaginary components of the scalar fields. This 
RGE then becomes, 

JTLT 

(47r)' = [2 {Xafbe + A'/.,e) ^ecf + Kbef^lef + iP ^ c)] 

+ Ae/6cHae/ + SA^^^j {^laf + H/ae) 

+ 2Tr {(mx - zm'^) [(vl {U^wj + V.Xf } + (X^t + ) X^Xf 

+ w^wjxf + I w„Wt + (X^ + X^) Xf } 
+W,*Uf V: + X^t I vjv: + Xi (X^t + Xf 
+ (6^c)]} 

+ 2Tr { (Wf (m;, + zm;^) { V^Uf + X^ W,* } 

+ (X^ + X^) (mx + tu^x) XfaXf + V, (mx + tu^x) WaUf 
+ Xi (mx + ^m:^) { W,Wt + (X^ + X^) Xf } 
+ V, (mx + ^m';,) X.^Vt 

+X^ (mx + ^m';,) { VJV: + Xj (X^t + Xf ) }) 
+ (6^c)} 

+ Tr {U^l (U^ + U^) + (U;^t + ) U;^, + VlV., + Wj^W, 
+\ {X^t + X^) + (X^t + x^t) XI, } } H.,,. 

+ Tr ju^mi + \J^M' + V^,V, + Wtw,, + 1 {X^txi + x^+X^, 

+ Tr ju^JU^ + U^^U^, + i {X^;X^ + Xf X^,}| (H,,,, + H:,,,) 

+ Tr ju^mj + + vt,V, + WtW,, + ^ {x^/X^ + Xf X,^,} | H.^'. 

+ Tr jujui + Uj, + i {Xftxi + xf Xj,}} (H,„, + H^O 

-3/ [C2(a) + C2(6)+C2(c)]H,b,, 

(13) 

where, as with (5) and (6), a sum over all gauge group factors is implied in the last term. The 
origin of the various terms in this long equation when compared with the two-component 
form of (7) should be clear. The terms involving the quartic couplings clearly originate in the 
first term of (7), the next set of terms involving the mass matrices for the Majorana spinors 
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in the second term, the terms with the summation over {a', b', c'} appear in both forms, and 
the last term involving the gauge couplings is clearly the corresponding term in (7). In the 
context of the MSSM, where we can without loss of generality choose for the index a in iiabc 
to always refer to a sfermion, the RGE can be somewhat simplified. This is because, when 
a is a sfermion index, U;^'^ and X^'^ are always zero, and the corresponding terms drop out 
from the equation. We remark that the first trace term not involving Majorana fermion 
masses in (13) is clearly symmetric under 6 c while the last two trace terms {i.e. those 
with c' as a dummy index) are simply the b ^ c equivalent of the previous two trace terms 
involving the dummy index b'. Thus the RGE for H^bc preserves the b ^ c symmetry of 
this coupling noted below Eq. (2). 

The RGE for the real scalar mass squared parameters in (1) reads [22], 



(4vr)' /3, 



1— loop 



ab 1 



(14) 



with 



Tiab =Tr { (y.Y,^ + Y,Yt) m^m} + (yIY, 
+Y„m}Y,m} + m/YtmyY;} , 

Ai = [C,{a) + C2{b)] ml, 



Kb 



Tr { Yl, Y, + Yt Y,, } m^, + Tr { Y^, Y, + Y^ Y,, } 



m 



ab' 



(15) 

(16) 
(17) 



where in (17) we have once again made the notation more explicit. Writing the complex 
fields as = i^ii±^ and comparing the coefficients of the bilinear terms in the real 
scalar fields (p, then gives. 



m 



B 



ab 



ab 

™ib + ^ab 



m 



ajbi 



(18a) 
(18b) 



In the i?-parity conserving MSSM, we never have non-zero entries in the Lagrangian corre- 
sponding to both ml and Bab for the same a, b. It suffices, therefore, to write the RGE for 
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just one of these combinations, which we take to be the first one. We then have, 
,2 d [mlf^ - Bab] 



[An] 



dt 



{2 [Xafbe + ^fabe) ™e/ " ^abef'^lf " {^efab + ^baef) ^ef] 



- 2 



2 Tr 



+ 



v^v: + w,w; + + x^) + x^ (x^t + xf )_ 

X (mx - imx) (mx + ^m^)} 
+ Tr { (X^ + X^) (mx + ^m';,) X^ (mx + ^m';,) } 
+Tr { (X^t + Xf ) (mx - ^m'^) Xf (mx - ^m;^)}^ 

3/ [C2{a) + C^m {ml -Bab) 
Tr {U^t (xji + U^) + (U^t + Uf ) U^, + VtV., 

+Wl,W. + 1x5 {K + + l + ) <'b 
- Tr {U^ (U^ + Vl) + (U^t + uf) \Jl + Vl,V. 

+WtW., + ix^t (x^ + X^) + ^ (X^t + ) X^} B,b 
+ Tr{u,^tu^ + Ufu2 +v;,V, 

+ ix^;X^ + ^X^X^j (m^,, - B,,,) 

+Tr (u^U^ + U^U^, + + ixf Xi I (m^,, - B^b'T 



(19) 



where, as before, a sum over gauge group factors in the term proportional to 6*2(5') is 
implied. Just as in the RGE of Eq. (13), we have the ordered the terms in (19) to make 
evident the correspondence with the terms in (14). The length of this equation is somewhat 
deceptive because many terms vanish in the case of the -R-parity conserving MSSM. In a 
similar manner to the RGE for Habc earlier, and as also seen in Paper I, when a and h 
are sfermion indices, U^'^ and X^'^ are zero. Likewise, when a and b label Higgs fields, 
^a,b and Wa^b are zero. The derivation of the MSSM RGEs is, therefore, considerably less 
cumbersome than it appears at first sight. 
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III. PARTICLE DECOUPLING 



In the previous section, we have obtained the RGEs for the Majorana fermion and scalar 
mass parameters, as well as the bilinear and trilinear scalar coupling parameters for a gauge 
field theory with spin-0 and spin-1/2 fields. While these RGEs are not applicable to a com- 
pletely general gauge field theory, as we have explained in Sec. II, they certainly apply to 
the MSSM with a conserved /2-parity quantum number. Since we did not assume supersym- 
metry in our derivation, it is straightforward to include SUSY-breaking threshold effects in 
essentially the same way as in Paper I. Specifically, we implement SUSY and Higgs particle 
thresholds as step functions in the evaluation of the /5-functions; i. e. we include the particle 
V in the effective theory only if the scale, Q, is larger than the mass of V. 

To implement particle decoupling using this procedure clearly requires a knowledge of the 
particle spectrum which, in many models, is obtained using the RGEs. Fortunately, because 
the results depend only logarithmically on the scale at which we decouple the particles, 
an approximate knowledge of the particle spectrum suffices in order to implement particle 
decoupling. Here, for reasons detailed in Sec. Ill of Paper I, we will decouple both higgsinos 
at Q = \fi\, the gauginos at the scale |Mj|, and the heavy Higgs bosons at the scale uih- 
The decoupling of sfermions, if mixing effects are negligible, is also straightforward. Since 
one of the main reasons for this analysis is to study flavour-physics in the squark sector [23] , 
we clearly must include mixing among the squarks. We will, therefore, defer the discussion 
of squark threshold corrections to a later point in the paper. 

It is evident that for any discussion of thresholds we need to know the mass parameters 
in the MSSM. In addition to the SSB scalar and gaugino bilinears given in (16) of Paper I, 
we have mass terms for the higgsinos as well as Higgs scalars from the superpotential. These 
higgsino terms are. 



resulting in two Majorana higgsino states of mass in the approximation that any gaugino- 
higgsino mixing can be neglected. 



scalar terms in the potential. The so-called D-term contributions to the Lagrangian density 




(20) 



The superpotential, /, also leads to Higgs scalar bilinears along with trilinear and quartic 



11 



also lead to quartic scalar couplings that, in the supersymmetric limit, are fixed by the gauge 
couplings. In the notation of Ref. [4], these scalar potential terms are given by, 

2 



1 ^ 

2 



-E 



df 



(21) 

s=s 



Superpotential interactions can result in scalar bilinear, trilinear and quartic terms. As an 
illustration, if we take the second term in (21) and choose to differentiate with respect to 
the up-type Higgs superfield, i.e. <S = Hu, we see that, 

C3- - U^jiku[i (f«)L {Q*im ^RmUL 



(/i%^");S«,/.° + h.c. 



(22) 



Note that we have inserted a tilde over the fi in writing these terms. This is to allow for the 
fact that the higgsino mass fi and the corresponding superpotential parameter fl in the scalar 
sector will, in general, evolve differently once SUSY-breaking threshold effects are included. 
This is the analogue of the corresponding situation for dimensionless parameters. As we 
saw in Paper I, gaugino- quark- squark couplings g'^ evolve differently from the corresponding 
gauge couplings (and in fact develop flavour- violating components), while higgsino-quark- 
squark coupling matrices f^^ evolve differently from the corresponding Yukawa coupling 
matrices iu^d once the renormalization scale Q is below the mass of the heaviest sparticle. 
We should emphasize that we can never get an RGE for the parameter fl discussed above. It 
enters via the SSB Higgs squared mass parameters in the combination ml ^ + or, as in 
(22) above, via combinations like (/i*f^")*^, the coefficients of "non-analytic" trilinear scalar 
interactions [24], denoted generically by Cj^ in Ref. [4]. It is only for these combinations 
of coefficients that enter the Lagrangian that we can (and do) obtain an RGE, not for the 
separate pieces. Above all thresholds, the RGEs for these combinations agree with the 
RGEs obtained from their component pieces. Notice that we have added a superscript $ 
to the 'Yukawa coupling constant" that enters the trilinear scalar interaction in (22). Here, 
where this term in the scalar potential originates in the derivative of the superpotential with 
respect to h^, we set $ = 

The reader may legitimately wonder why we include a superscript on the coupling f,'^" 
in the trilinear scalar coupling term but write the quartic scalar coupling constant as the 
square of the usual quark Yukawa coupling even below the scale of SUSY breaking. The 
reason is that we are going to ignore the difference in the renormalization of the quartic 
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scalar couplings, which is also why we did not derive the corresponding RGEs in Paper I. 
These couplings are less important from a phenomenological perspective (even though some 
of them enter the squark and slepton mass matrices). It is possible that quartic squark cou- 
plings that arise from the first term in (21) may also develop a non-trivial flavour-structure. 
These couplings will not affect the two-body (flavour-violating) decays of squarks except 
at the loop-level, leading us to believe that threshold corrections to these couplings can be 
sensibly neglected. In the following, we have set all quartic scalar couplings, irrespective 
of whether they originate in the superpotential or in the D-teim, to their supersymmetric 
values, i.e. equal to the "square" of the usual Yukawa coupling or of the corresponding 
gauge coupling, but have retained threshold effects in the trilinear couplings that cause 
phenomenologically important mixing between left- and right-squarks.^ Since we do not, 
therefore, need independent RGEs for these quartic couplings, we have not listed here equa- 
tions to convert the A and A' in (2) to the corresponding A in (1), equivalent to (12) and 
(18). Of course, the inverse of these equations is needed in order to write the right hand 
sides of (13) and (19), but this is considerably simpler. 

The decoupling of fields whose mixing can be ignored is straightforward: below their 
mass scale, we simply remove all contributions from these fields when evaluating any RGE. 
This covers: gauginos; higgsinos, because they are approximately degenerate [14]; and in 
many models, also sleptons. 

Next, we turn to the decoupling of the spin-zero particles in the Higgs sector where we 
cannot disregard mixing effects, and then discuss how we treat decoupling in the squark 
sector where evaluation of flavour effects is one of our important goals. Before proceeding, 
we mention one more complication that arises because fi does not necessarily equal fl. This 
makes the determination of fi more involved since the electroweak symmetry breaking con- 
ditions only depend on its bosonic cousin fl; we will defer the discussion of how fi and fl are 
determined to the next section. 

^ The quartic terms in the potential yield the so-called D-term contributions to squark masses, and also the 
term that makes the squark mass the same as the quark mass in the supersymmetric limit. In most cases 
these are both subdominant contributions to squark masses, giving futher credence to our approximation. 
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A. Higgs Boson Decoupling 



As discussed in Paper I, the implementation of step-function decoupling requires us to 
write the Lagrangian density with the fields in their (approximate) mass eigenstate basis. 
This led us to rewrite the interactions of the MSSM Higgs fields in terms of the spin-zero 
fields (h, H^) defined [14] in terms of the MSSM fields of the scalar Higgs sector by, 



(23) 
























hi) 


\h'X 



(24) 



where c = cos P and s = sin p. For ttih ^ Mz, the physical Higgs bosons h, H, A and if^ 
are then approximately given by [14],^ 

^ ~ V2 



n 



-H+iA 
V2 



Here and are the would-be-Goldstone bosons that get dynamically rearranged to 
become the longitudinal components of the massive gauge bosons by the Higgs mechanism. 
If CP is violated in the Higgs boson sector, the neutral spin-zero states H and A will further 
mix with one another. When we eliminate the charged and neutral components of the 
scalar doublets and Hd and rewrite the Lagrangian density in terms of (h,7i, H^), new 
operator structures, and concomitantly new couplings, appear in the theory. For instance, 
the couplings of h to fermions are given by the SM Yukawa coupling matrices A^ ,;;, and 
after decoupling of the heavy Higgs scalars, the matrices disappear altogether (though, 
depending on the spectrum, their tilde cousins may remain). 

For Q > ntH where all Higgs fields are active in the RGBs, the rotation to the Higgs field 
"mass basis" makes no difference as it is just a field redefinition. For Q < mn, however, 
this rotation is crucial since it determines the particular combinations of fields that decouple 
from the RGEs. In this case, only specific combinations of the parameters in the original 



^ If rriH is close to Mz , all the bosons in the Higgs sector have masses close to Mz and, as noted in Paper I, 
the resulting threshold corrections are small. 



14 



theory remain upon decoupling, while other combinations generally become irrelevant. For 
example, when the two trilinear terms which couple ul and ur to the Higgs bosons are 
rotated into the (h,?i) basis, we see that 



= U^Rki^DklULl {Sh + CH) - M^jfc ULl (ch - SH) . 

Upon decoupling the field 7Y, this term becomes 



(25) 



s(a^),,-c(/i*e)^^ uuh. (26) 



Although in the complete MSSM case we can sensibly talk about the evolution of the 
constituent pieces (a^)^^ and fi separately, for Q < mn, it is no longer possible to write 
RGEs for both and /i*f^". We must instead talk only about the single combination 
s{a^)ki — c (/i*f^")^^ that remains in the effective theory below Q = mn- In the same 
vein, we should mention that the quartic scalar operators with A'^^^^-type couplings in (2) 
are absent above all SUSY thresholds in the i?-parity conserving MSSM, since in this case 
the number of daggered and undaggered fields is always the same (remember that these 
originate in the absolute square of a quadratic operator). We will see below that these 
A'-type couplings arise because of the way we define specific linear combinations of fields to 
take into account threshold effects. 



B. Squark Decoupling 

Squarks of different flavours and types (L or R) can mix as long as they have the same 
electric charge. In the up-squark sector, squark mass eigenstates are, therefore, combinations 
of ul,r, cl^r and ii^R squarks, and likewise in the down-squark sector. It is these mass 
eigenstate flelds that need to be decoupled for the evaluation of threshold effects. 

The reader may have noticed that in Paper I, as well as in this paper, we have maintained 
manifest electroweak gauge invariance throughout, in that we decouple entire multiplets to- 
gether. In keeping with this, we have specifled the boundary conditions for the (3 x 3) 
Yukawa coupling matrices at the scale Q = rrit, rather than a.t Q = Mz where the top 
quark would have already been decoupled.'^ To avoid any complications that may arise 

^ Specifying the boundary eonditions at Q — irit is really only necessary if we insist on decoupling every 
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from splitting a multiplet, for the purposes of evaluating the squark thresholds, and for 
this alone, we will assume that left-right squark mixing effects which can arise only from 
SU (2)l X U{1)y breaking are not important.^ This is clearly a sensible approximation when 
squark SSB parameters are larger than the weak scale. Remembering that the calculation 
is only logarithmically sensitive to the actual location of the threshold, we can see that this 
approximation breaks down significantly only when the off-diagonal entries in the squark 
mass matrices cause large cancellations in the calculation of the lighter squark eigenvalue, 
leading to a physical mass squared much smaller than either diagonal entry. For patho- 
logical parameter values where such a cancellation is operative and ql — Qr mixing indeed 
causes one of the squarks to be much lighter than all other squarks of the same charge, our 
approximation will not apply. Hereafter, we will assume that any hierarchy in the squark 
mass spectrum has its origin in the values of squark SSB parameters, and is not the result of 
an accidental cancellation between the diagonal and off-diagonal parameters in the squark- 
mass matrix. We will, therefore, neglect left-right mixing, for, and only for, the purpose of 
determining the locations of the squark thresholds. 

Squark thresholds are then determined only by the SSB mass squared matrices for the 
squarks. As with all other particles, we decouple any squark below the scale equal to its 
mass. Since (for the purpose of locating the squark mass thresholds) we are working in 
the approximation that SU{2) x f/(l)y breaking terms are neglected, the various left- and 
right-squark mass thresholds are determined by the eigenvalues of the corresponding SSB 
mass squared matrices. The associated technicalities are best explained by describing the 
step-wise procedure that we use. 

1. We evolve the gauge and Yukawa couplings to the high scale (usually taken to be 
Mqut^) where the boundary conditions for the SSB parameters are specified. Toward 
this end, we first eliminate the quark masses in favour of the diagonal SM Yukawa 
coupling matrices in the quark mass basis at Q = Mz {Q = rrit for the top quark), and 
rotate to a current basis (related to the quark mass basis by Eq. (38) of Paper I) at Q = 

particle at the scale equal to its mass. Such a procedure is advantageous in many situations, but not 

compulsory. Indeed, as we will see in the next section, while we decouple all sparticlcs and heavy Higgs 

bosons at their mass scale, we decouple the top quark at Q = Mz, below which the effective theory is an 

SU{3)c X U{l)em gauge theory. 
^ To be specific, we will focus our discussion on squark thresholds, but exactly the same procedure applies 

for slepton thresholds. 
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rrit.^ We morph the SM coupling matrices Xu,d to the MSSM Yukawa couphng matrices 
iu,d at Q = iTiH, and include SUSY radiative corrections [25] as given by ISA JET [26] 
{i.e without flavour mixing among squarks) at Q = MgusY, before continuing to evolve 
to Mqut- 

2. All the dimensionless couplings and MSSM SSB parameters can now be evolved back 
to the low scale with Higgs bosons, higgsino and gaugino thresholds implemented as 
in Paper I. We must, however, be careful in implementing the squark thresholds. At 
each step in the integration, we flnd the eigenvalues of the SSB squark mass squared 
matrices, and carry on the evolution in the original current basis if the scale is 
larger than the highest of these eigenvalues, and none of the squarks are decoupled. 

3. At some value of the scale, Qq, just crosses the largest eigenvalue of one of the 
SSB squark matrices, evaluated at the scale Qq. We decouple that particular 
squark at this scale, but retain all other squarks in the evolution. The up and down 
type squarks in the same doublet decouple together of course, so that SU (2) is never 
broken by this procedure. To decouple the squark, we rotate to the basis (denoted 
here by the superscript M) where the corresponding squark SSB matrix is diagonal at 
Q = Qo using the appropriate one of, 

(27a) 
(27b) 

RJ^ , or one of (27c) 

(27d) 

cr = Re^R when decoupling a slepton. (27e) 

The unitary rotation matrices, R,, are chosen to diagonalize the Hermitian SSB squark 
(or slepton) mass matrices, for example: 

(R^m^R^)^^ = (u^^ ^ in^t'S., • (28) 




5 



In our iterative procedure to solve the RGEs described below, it is important to remember that we retain 
the top quark all the way to Q = Mz so that electroweak gauge invariance is preserved. 
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Decoupling this squark is now straightforward — as with all sparticles, we simply 
introduce a 6q^, for the decoupled squark g^'^ (where k is the mass basis index, and 
Q = Qli or dfi) into the RGEs, since we are in the "mass basis" for this squark at 
least within the approximation that we have discussed.^ Just for clarification, in any 
other basis these would be matrices. For instance, the left-handed squark theta, 
6q^, can be written in the original current basis as, 

(0q),, = (Rq&'^^^R^q),, = eQ^{RQ),MK)i, . (29) 

We have checked that this same 0q matrix works for all terms involving doublet 
squark thresholds, and likewise for the other squarks and sleptons. We always write 
the RGEs in a current basis where the SSB squark mass matrices are diagonal (see the 
item immediately following this). This is why we have 6's rather than the matrices 
appearing. We rotate back to our original current basis when we present numerical 
results. 

4. To evolve to lower values of Q, we must be in some current basis. We thus also rotate 
the corresponding quarks (or leptons for slepton decoupling) by the same amount as 
the decoupled squark (slepton), and continue the evolution in this new current basis. 
The "mass" and the "eigenvector" of the decoupled squark is frozen at its decoupled 
value at Q = Qq. As a technical aside, we note that the same RGE (with (3 x 3) mass 
squared matrices) can be used to continue the evolution in this new current basis. 
The Oqf,'s that we introduce into the RGEs ensure that the decoupled squark does 
not contribute to the evolution of the remaining squarks. In other words, in the new 
mass basis the decoupled squark never contributes to the evolution of the SSB matrix 
elements for the (2 x 2) sub-matrix in the sub-space orthogonal to the decoupled squark 
state.'' We use the RGEs to continue the evolution to lower scales, now obtaining the 
eigenvalues of the (2 x 2) sub-matrix in the orthogonal sub-space just mentioned, and 

^ From now on, we will no longer write the superscript M for the squark mass eigenstates, but it will 
presumably be clear from the context whether we are referring to these, or to the states in the original 
basis. 

^ Rather than rc-code new RGEs where we retain only the two active squarks for the subsequent evolution, 
we continue the evolution with all the squarks. Below the scale Qo, just the elements of the (2 x 2) 
sub-matrix in the vector space orthogonal to the decoupled squark are physical. The remaining elements, 
though calculated are never used, and their value does not affect the calculation. 
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once again decoupling the squark when the larger of the two eigenvalues just crosses 
. The corresponding eigenvector of this sub-matrix gives us the second squark state 
with the same gauge quantum numbers to be decoupled (remember that we are in the 
new current basis). As before we freeze its mass and wave function upon decoupling. 
Notice that, by construction, the two decoupled states are orthogonal to one another 
and to the single remaining active state, as they should be. The single remaining state 
may be decoupled in a straightforward way. Since we know the relation between the 
new current basis and our original current basis, it is straightforward to rotate the 
SSB parameters back to the latter basis. 

5. We have described squark decoupling for any one set of squarks, i.e. doublet squarks, 
the up-type singlet squarks or the down-type singlet squarks, but it should be clear 
that the procedure that we have described works for all the squarks as well as for 
sleptons. 

6. We must, likewise, decouple squark mass eigenstates obtained using our procedure 
during the evolution of all other parameters. However, for the scalar trilinear coupling 
parameters, a, there is the added complication that we must decouple the heavy Higgs 
bosons at the scale Q = rriH. In this case particular combinations of parameters that 
include need to be evolved, as we have already discussed. 

IV. APPLICATION TO THE MSSM 

We now apply the general results obtained in Sec. II to derive the RGEs for the dimen- 
sionful couplings of the i?-parity conserving MSSM. As in Paper I, our strategy will be to 
read off the couplings in (2) from the MSSM couplings detailed in Eq. (16)-(22) of Paper I, 
and from (20)-(22) of this paper, and substituting these into the general RGEs obtained in 
Sec. II. The quartic scalar couplings not listed here may be derived from (21) above, and 
are also found in Ref. [4]. Once the couplings are all correctly identified, the derivation of 
the RGEs using (5), (6), (13) and (19) is tedious but straightforward. 

We proceed by outlining how to go about deriving the RGEs for each group of dimen- 
sionful couplings: gaugino and higgsino mass parameters, trilinear couplings, and scalar soft 
mass parameters and their bilinear cousins. Our purpose is to guide the reader interested 
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in deriving the RGEs as to where the various terms come from, and to point out potential 
pitfalls that may be encountered on the way. The complete set of RGEs for the dimensionful 
parameters is listed in Appendix B. 

A. Gaugino and Higgsino Mass Parameters 

The application to the MSSM of the general RGEs in (5) and (6) is straightforward. The 
matrices m^'' are constructed from the m['\^, the SSB gaugino mass terms, and /i which 
appears via the higgsino mass terms in (20). The parameter fl that appears in the scalar 
Higgs sector should not be confused with the corresponding fermion mass parameter /i. 

As previously noted, the trace terms in (5) and (6) vanish so that we need only sum 
over the scalar index b and multiply out the terms. When the higgsinos are rotated as in 
Eq. (14) of Paper 1, m^^ becomes diagonal and the derivation simplifies further. We find, 
for example, that the RGE for M2 is 




+ 2sc (-0, + Oh) Qi, [g'''li*t- + (/'')>(^'")*] - l2e^M,gl , 
and draw the reader's attention to the following points: 

1. The first term, which is M2 multiplied by a number of gaugino coupling terms, arises 
from the terms in (5) that have (mx ± zm^) on the extreme left or on the extreme 
right. When in this position the m^^ is connected to an external gaugino, and since 
m^'' is diagonal it contributes a gaugino mass term. This term includes a sum over 
scalar fields and corresponding fermion fields (that enter the RGE via 1-loop Feynman 
diagrams), with concomitant couplings such as g*^ and g^^ of fields that couple to the 
wino. This sum over all scalar fields includes all active sfermion flavours at the scale Q, 
i.e. sfermions with masses larger than Q are decoupled. As in Paper I, the decoupling 
of any particle is introduced into the RGEs via the 6*^ = 1 if Q is larger than the mass 
of "P, with Op = otherwise. Thus, for instance, in the first term 6'^^ (g'^)fci(g'^)|^, the 
summation over k includes all active left-type squarks at the scale Q (as emphasized 
in the last section, this RGE is written in the basis where the corresponding squark 
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SSB matrix is diagonal) , while the / sums over the quarks (that are all always assumed 
to be present in the effective theory we use to compute the RGE*) which is why we do 
not include the corresponding ^'s. Below the highest sfermion threshold the sfermion 
gaugino coupling terms are, therefore, truncated traces. 

2. In the RGEs, we write 6'p's for, and only for, particles V that enter via the internal lines 
of loop diagrams that contribute in the RGEs. Sometimes, this internal particle is the 
same as the particle on the external leg of the corresponding diagram, as illustrated 
by the appearance of 6^r in the first and third terms of (30). Furthermore, since 
at any scale Q, we retain only those (non-SM) particles with masses larger than Q 
in the effective Lagrangian that we use to derive the RGEs, we stop evolving the 
coefficient of any (composite) field operator — these coefficients are just the Lagrangian 
parameters — when the scale Q falls below the mass of any of the fields that enter 
the operator. Thus, in the case of M2{Q), we freeze its evolution below the scale Qo 
where Qo = ^2(Qo)- 

3. The reader may be struck by the second term on the right-hand-side, proportional to fi, 
whose appearance seems odd at first sight. It originates in the terms in (5) where m^^ 
is sandwiched between two X^'^ matrices. Since X^'^ connects gauginos to higgsinos, 
and the external fields are gauginos, the 111^ is necessarily a higgsino mass term. 
Notice that this term is a threshold effect: it does not vanish only ii \fi\ < Q < mn- 
The light and heavy Higgs boson doublets in (23) and (24), respectively, make equal 
but opposite //-dependent contributions to the evolution of M2 which indeed cancel 
above all thresholds. The appearance of /i is a general feature of the electroweak 
gaugino mass RGEs — listed in full in Appendix B, (B2)-(B5). For much the same 
reasons, for appropriate mass ordering, the RGE for the complex parameter /i, (Bl), 
develops a dependence on the electroweak gaugino mass parameters, m['^ and Mg'''. 

4. Despite the appearance of the complex //-dependent terms just discussed in the RGEs 
for the gaugino mass parameters, the reality of Mj and M'- under renormalization 
group evolution is preserved, as it must. 

In this paper, we always assume that Q > Mz, and retain all SM particles in the effective theory, but 
write an explicit Oh only for the light scalar doublet that includes the would-be-Goldstone bosons. 
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5. Finally, notice that above all thresholds, where the SU{2) gaugino coupling matrices 
g* reduce to g2 times the unit matrix in the flavour space (see Paper I) and all 6i = 1, 
we recover the MSSM result upon summing over all flavours [18]. 



B. Trilinear Couplings 

For the i?-parity conserving MSSM, the trilinear scalar couplings Habc involve couplings 
between one of the Higgs boson flelds and two sfermion flelds. In order to facilitate de- 
coupling of the Higgs scalars, we saw that it is necessary to rotate these Higgs flelds to 
their (approximate) mass basis given by (23) and (24), so that a, b and c run over the com- 
plex fields, |h, 7Y, G"*", iy'^, Mij, Jlj, e^i, i/Li, Mffi, c/^^ and i runs over all three flavours. 
Unlike the charged would-be-Goldstone flelds that appear explicitly, the neutral would- 
be-Goldstone boson is contained in the complex fleld h. When fully expanded out in 
flavour space, Habc is a (25 x 25 x 25) array, whose entries are mostly all zero, and with 
many non-zero entries related by the SU{2) gauge symmetry of the interactions, that we 
can easily read off. For instance, from (26) we immediately see that, 



H 



si^lh-c{ji*i!:-)l . (31) 



As we mentioned below (13), we can without loss of generality in the RGBs, always choose 
the external index a to be a sfermion, since then several terms drop out in the derivation of 
the corresponding RGE for the trilinear parameter. Of course, we also need to work out the 
A and A' matrices. These are (25 x 25 x 25 x 25) arrays in scalar fleld space, again with 
mostly zero entries and with most non-zero entries related by the gauge invariance of the 
interactions. As noted previously, these can be worked out from (21). The extraction of 
the A terms is entirely straightforward, but care must be taken with the sfermion flavours 
and squark colours, which are implicitly summed over in (21). As alluded to at the end of 
Sec. HI A, we note here that A' = when the MSSM Lagrangian is written in terms of hu 
and hd Higgs flelds. These couplings arise only when the Higgs boson doublets are rotated 
to their mass basis using (23) and (24) since only the conjugate flelds h^* and h^* appear 
in the linear combinations with the unstarred flelds and h^, respectively. In order to get 
only one daggered scalar fleld, lepton and baryon number conservation for the dimension 
four operators imply that we must have Higgs flelds in the interactions. Moreover, we must 
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have both an up-type and a down-type Higgs field before the Higgs field rotation to the mass 
basis, and that the down-type Higgs fields must he hy while the up-type fields must be 
^(+,0) r^Yiis only occurs in quartics in (21) that derive from differentiating the superpotential 
with respect to the superfields ul ot d^, and so have a and a dj^ in the interaction. As 
a result, there are very few no n- vanishing A'-type quartic interactions. Indeed, these are 
completely given by, 

^ 3 iO:) ,1 i^H^ - G^n) u^^dni + h.c. (32) 

In deriving the RGEs, we must remember to include the contributions from the would- 
be-Goldstone fields in the sum over complex scalar fields. The neutral Goldstone boson 
is automatically present in the field h along with the (almost) SM-like Higgs boson, but 
the charged fields must explicitly be included in the sum. For example, the ROE for 
'^umXuLj includes a term Kej^h^uLj^uRi,ej and we must sum e and / over all scalars in the 
theory for which Aej,h,«i,H,a^^,ej ^ 0. 

Using Eq. (13) to derive the RGE for the operator in (31), writing all internal threshold 
0p's explicitly, we obtain: 
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(47r 



2d|^s(au)ij-c(/i*f^' 
' dt 



kj 



+eHSC [c(a„),fc + s 



^9 



12 



-2(^ + ^)M,+6(y,,(f„) 



g 4, + 2[(f.)t(f„)]^,^. 

[s(a„)fc,-c(/i%^")J 



+ 



X 



s(a„)fc,-c(/i*C")^^. 



Ohsc 



Sal 



5,, + 2[(f„)(f„)t] 



+ [(f,)(f,)t]J X [c(a„),, + s (rt") 



kj 



.(a,),fc + c /i*f^ [(f,)t(f„)] 

V / jfcJ 
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^^.(M3-zM^)(g«)*,(fJ,,(gr 



- 6e^e-,s (M, - iM',) {g^-n^'^)UiT)k, + fBO-HCf^^g'^' (4(f«),,(g--)^^. - (g'^)*,(fr 

+ 6%%c/i*^'^^(g«)*,(C«),, - 4^^c/i*(f«),,(f])H(fr 



3s2(ft),,(fJ,, + {3(f])fcKfd);fc + (fl)fc/(fe);A:} 
+ Ohsc \?>{fl)ki%)ik - {'i{il)ki{U)ik + (fl)fcKfe) 
+ 



s(a„),,-c(/i*e), 



X 



s(aji, - c (/i*f„^") 
(a„),, + s(/i*e")^^. 



+ en. 



s(a„),, -c(/i%^"),. 



(33) 



The first set of terms before tlie appearance of (Mi — iM[) arises from tlie A and A' terms 
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in (13). All A terms have a non-zero contribution but since most A' entries are zero, only 
Kbef'ti*cef contributes. 

The second set of entries, containing gaugino mass terms, originate in the non-vanishing 
traces in (13) where [nix — i^'x) is on the extreme left. For our case, $a = u/jj, so that 
Vq has a quark/gaugino as the first /second index, while has the higgsino/quark as the 
first/second index. Remembering that taking the transpose, or the dagger, flip the order 
of these indices, it is quite straightforward to see that the first set of trace terms where 
(mjf — i'm'x) is located on the extreme left necessarily give contributions proportional to 
SSB gaugino masses. For this, we must keep in mind that the matrices X^'^ connect gauginos 
to higgsinos, and vice- versa. A similar analysis of the next set of trace terms in (13), with 
(mx + iT^'x) enclosed by other fermion matrices, shows that the altered location of the 
Majorana fermion mass matrix now results in the appearance of the "higgsino mass", /i* 
(without any tilde). 

The remaining terms, except for the very last one which obviously arises from C2{S)liabc, 
are from the traces with primed scalar indices. The important thing to note is that the trace 
in (13) denotes a sum over fermion types, which sometimes, but not necessarily, becomes a 
trace over fermion flavours, since now a scalar may carry a flavour index. In the present case, 
the trace over the product of two U^'^ matrices does result in a trace over fermion flavours, 
since then the scalar index on the U^'^ necessarily corresponds to a Higgs field which does 
not carry any flavour. We refer the interested reader to the discussion below Eqs. (36) and 
(37) of Paper I where a completely analogous situation is discussed in more detail. 

We have illustrated the derivation of the RGE for the trilinear scalar coupling of squarks 
to the lighter of the two Higgs doublets. The RGE for the corresponding coupling, 
c{au)ij + s (/i*f^"')^_^. , to the heavier doublet can be obtained in the same manner. By 
taking linear combinations of these RGEs, we can obtain the separate RGEs for a„ and 
/i*f^".^ Since the coupling /i*f^" always occurs as a product, it is not possible to obtain the 
RGEs for the individual factors. Of course, above all thresholds, we must have /i = /x and 



^ The s and c can be taken out of the derivatives on the left-hand-side as noted in Paper I. For Q > mn, 
when both doublets are in the theory, the rotation is irrelevant. Otherwise, the input value of tan/3 
corresponds to the ratio of VEVs at the scale Q = rriH when the two doublet model reduces to the one- 
doublet model. We may think about this as evolving the couplings a„ and f^" from the high scale down 
to the scale Q — niH, at which we must do the Higgs rotations to reduce to the one-doublet model. 
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fhu — have checked that with these replacements, our RGEs reduce to the MSSM 

RGEs [18] if we put all 6i = 1 and take care to sum over all internal flavours. 

C. Soft Masses 

Finally, we turn to the RGEs for the SSB mass parameters and their bilinear cousins B; 
we have almost all the matrices necessary to use (19) to find the RGEs. The majority of the 
missing and B terms appear in the SSB Lagrangian and so the required matrices can be 
written down directly. Note that when we write the Lagrangian in the rotated Higgs basis, 
there are no 13 terms so that Bg/ = 0. In this case, several terms drop out of the RGE in 
(19). 

1. Higgs Mass Terms 

Using (19) to derive the RGE for the coefficient of the h''"h term in the Lagrangian we 
find that 
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's'\~g'H^ + c'\g"'A^ 
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's'\~g'H^ + c^\~g'^A^' 



+ [Ml + M'i) 
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{-49-^9 ^scfi*~g''^-g'^'^ (Mi + iM[) - 129-^9^scfi*~g^-~g^^ {M^ + zM^)} 
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i {-4^^^^sc/i(^'^")*(^"^'')* (Ml - zM[) - 129-,9^scf,{~g'^ri~g''r (M, - tM',)} 
9u (^ + ^) [s'Ml^ +c'Ml^-scib + b* 
0^ { [eCC] + [6f,*fJ + 2f:fJ] + 9^9-^ 

[s'M'^^+c'Ml^-sc{b + b 
Ohsc \ [QrX]^^ - [6f,*fJ + 2f:fJ],,, + 9,9^ 



s'\~g'^A^ + c'\~g"'A^ 



}L{Mi-M|J-i(c^-.^) (6 + 6* 



(34) 



where M^^ = (m?,^ + Iftf) and M^^ = {ml^ + 

All terms up to the trilinear scalar couplings derive from the single term, Aa/bei^lf, in 
(19). All other quartic terms are zero, either because A' vanishes when a = 6, or because 
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Be/ = 0. Since this operator has no flavour indices, all flavours are internal and, therefore, 
summed over. Extra care should be applied when dealing with sums over squarks since there 
are additional factors of 3 arising from a sum over colours. Once the entries of the A matrix 
have been worked out, this contribution should not pose any special difficulty. The so-called 
S'-term is also contained in the terms proportional to g'"^ from this contribution. 

The terms containing the trilinear scalar couplings obviously come from the 
(H*^j + H/ae) Heftj term in (19). Since we have chosen the first index of ^abc to always 
be a sfermion, 'Raef = for a = h. For the contributions from the squark-Higgs scalar 
trilinear interactions the factor 2 changes to a factor 6 due to a colour sum. 

Next, we turn to the contributions from the Majorana fermion mass terms from the m^^ 
matrices in (19). When the two factors of (mx ± zm^) are next to one another, it should 
be clear that we should obtain either a or an (Mf 2 + M'^^ term. On the other hand, 
when the factors of (mx ± ^m'j,^^) are separated by an X^'^-matrix, we obtain a product of 

1 2 

gaugino and higgsino mass parameters because X^' only connects gauginos with higgsinos. 
The trace, of course, is a sum over all gauginos and higgsinos in the theory. 

The terms that derive from the next term, C2(5')m^{,, should be obvious. Following these, 
we have the set of traces with primed scalar indices, a' and h' . As discussed earlier, many of 
these terms are zero when a and h are Higgs scalars. We are left with traces over only Ui'^ 
and X^'^. The trace over the U^'^ matrices leads to a trace over matter fermion fiavours, of 
which the quark traces acquire an additional colour factor of 3. 

The Higgs scalar mass term for which we just obtained the RGE written above is the 
only combination which remains in the effective theory for Q < nin- The complete set of 
mass terms in the rotated Higgs basis are, 

C3-[s^ {ml^ + Ifif) + {ml^ + \fif) -sc{b + b*)] h+h 

- [c' [m\^ + l/iH + [m\ + + (6 + \f)\ U^U ^^^^ 

- {sc {m\^ + - sc {mjj^ + \fif) - c% + s%*] h^T^ 

- [sc {m]j^ + - sc {m]j^ + + s% - c%*] H^h , 

and the RGBs for these coefficients can be obtained in a similar manner. We can then obtain 
the separate RGBs for the real parameters {fn'jj^ + {itt-'h^ + and the complex 

parameter b, valid for Q > nin, by taking appropriate linear combinations. 

It is important to note that since the terms rrijj^ and only appear in the Lagrangian 
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in the combination (tti^^ + \ flf), we cannot derive an RGE for them separately. Of course, 
above all thresholds, supersymmetry requires /i = /i, so that we can use the RGE for /i, 
(Bl), to extract the RGE for the soft mass parameters m'jj^ and m'jj^, but these will cease 
to be valid once any one particle is decoupled from the theory. We have checked that these 
RGEs reduce to the standard ones [18] once all 6i are set equal to unity, and that our RGE 
for reduces to the RGE for the SM Higgs boson mass parameter [27] if we set the quartic 
scalar coupling in the SM to be the appropriate combination of gauge couplings. 

2. Sfermion Mass Terms 

The RGEs for sfermion soft mass terms are somewhat simpler on account of the fact that 
they have no Higgs fields in their operator in the Lagrangian. Otherwise, the derivation 
follows in a similar fashion to the Higgs mass terms, with the obvious differences in the 
terms in (19) that contribute to the RGE. There is still only one A term (because again 
the A' and 13 are zero), and this contributes a large proportion of the RGE including terms 
with, and terms without, traces over sfermion mass matrices. These terms include the S- 
term as for the Higgs mass RGE above. When e and / are both Higgs fields, this gives a 
6-parameter dependence in the RGEs for the sfermions listed in Appendix B for Q < mn- 
Since the external fields are sfermions, both trilinear terms in (19) now contribute, and in 
the trace terms, only V, and W, terms are non-zero. These terms yield the contributions 
that depend on the square of the gaugino mass parameters, and also on Above all 

thresholds, these latter contributions cancel with jl^ contributions arising from Higgs boson 
loops. 

There is one point about the RGEs for the sfermion SSB mass parameters that we ought 
to draw attention to. For Q < mn, where we only have the light Higgs doublet in the 
theory, we would expect that the couplings can occur only in the combination sf^. A 
look at the term just before the trilinear coupling terms in (B25) shows, however, that this 
is not the case. (There are analogous terms in the RGEs for mfj^,, and also for m| ^.) The 
"product of Yukawa couplings" that appears in these terms (without accompanying or 

factors) is really a Aikij-tjpe coupling for the quartic interaction of squarks which, as we 
have already explained, we have approximated by its supersymmetric limit and set equal to 
the product of the corresponding elements of the Yukawa coupling matrix, each frozen at 
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the scale Q = m^,; e.g. Aq^^^^-^^q^ ~ (: 



f'u)ifc (fn)^ ill the first of such terms in (B25). 



We have checked that, except for terms involving couplings and mass parameters of Higgs 
boson fields, the RGEs that we obtain agree with those in Ref. [17]. Indeed the RGEs in 
Ref. [17] do reduce to the MSSM RGEs if we set all the 6's to be one. However, since the 
RGEs of Ref. [17] appear to have been written without any rotation of the Higgs boson fields, 
we found it impossible to compare contributions involving thresholds for Higgs boson fields. 
For these same reasons, we are unable to see how their RGEs for Yukawa couplings reduce 
to the corresponding RGEs for the SM [14] when all new particles are decoupled. Likewise, 
we are not able to obtain the RGE for the SM Higgs boson mass parameter using the RGEs 
for the Higgs scalar SSB parameters as given in Ref. [17]. 

V. SOLUTIONS TO THE RGES AND FLAVOUR- VIOLATING SSB PARAME- 
TERS 

While the quark mass basis, where both up and down quark Yukawa matrices are diagonal 
(at a chosen scale), may be the most physical basis to work in, as we saw in Paper I, it is 
more convenient to work in a current basis where either the up or the down, but not both, 
Yukawa matrices are diagonal. In this case, we rotate the entire quark doublet (in the flavour 
space) thereby preserving the SU{2) gauge symmetry. To similarly preserve the underlying 
supersymmetry, we should also rotate the squark multiplets the same way that we rotate 
the quarks; i.e. we should rotate the quark (and also lepton) superfields. Of course, just 
as the Yukawa coupling matrices transform under this change of basis (see (38) and (39) of 
Paper I to set the notation), the SSB mass and a-parameter matrices in a arbitrary basis are 
related to the corresponding matrices in the basis where the up (or the down) quark Yukawa 
coupling matrix is diagonal aX Q = rrit (the matrices in this special basis are denoted with 
the superscript M) according to. 



where we set q = u{d) in ^ in the basis where up (down) type quark Yukawa coupling 
matrices are diagonal at Q = m^. The matrices Y ^^^{u^d), defined via (38) of Paper I, 




V«(«,rf)(a*yW[(g), 
V^(g)K)*V[(g), 
Vij(M,rf)(m^_^)^V{j(M,d) , 



(36) 
(37) 
(38) 
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are the linear transformations that connect any current basis to the quark mass basis. The 
SSB matrices and (ml)'^^ in the basis where up quark Yukawa couphngs are diagonal 
are related to the corresponding matrices in which the down quark Yukawa couplings are 
diagonal by, 



Here u and d in the parenthesis denote whether the up or down type Yukawa coupling 
matrices, respectively, are diagonal at Q = rrit, after the complete transformation (38) of 
Paper I; i.e. when we go to the basis where the up-type Yukawa matrix is diagonal, we also 
transform the singlet down sector by the matrix Vi?((i), and likewise if we transform to the 
basis where the down-type Yukawa coupling matrix is diagonal. The matrix K that appears 
in (39) and (40) is the Kobayashi-Maskawa matrix [28] given by. 



There is no a priori preference for the one or the other basis. However, if in a process 
there are only up, or only down, type quarks in the initial and final states, it makes sense to 
use the mass basis for these quarks. In the next section, we will apply the considerations of 
this paper to the rate for the decay of the lightest up-type squark to the final state cZi. For 
this reason, other than where explicitly stated, we will work in a basis where the up-type 
quark Yukawa couplings are diagonal aX Q = rrit, i.e. we will take y^L,Riu) = 11. We then 
have Vi((i) = K. Finally, to fully specify a current basis we must also choose Vj:}((i) — the 
matrix for the rotation of the right-handed down- type quarks. Except when we illustrate 
results for extensions of the mSUGRA model, we also set Y^^d) = 11. 

Since the boundary conditions for the dimensionless couplings of SM particles are specified 
at the weak scale, while those for the dimensionful parameters of the MSSM are usually given 
at a high scale, we solve the RGEs using an iterative procedure. This is, of course, quite 
standard, but the incorporation of threshold effects entails some new complications. We 

Within the SM, physics is independent of our choice of the matrices Vi?(u) and Vij((i). and depends on 
Vi(u,(i) only through the KM matrix. This is also the case for the MSSM with mSUGRA boundary 
conditions, but more generally, physics also depends on the right-handed quark rotation matrices. 




(39) 
(40) 
(41) 






(42) 
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begin by evolving the measured gauge and Yukawa couplings of quarks and leptons, specified 
in a chosen current basis, from Q = Mz to the high scale (usually Q = Mgut) where the SSB 
parameters of the MSSM are specified. We can now evolve all the MSSM parameters down 
to the weak scale, decoupling the particles one-by-one as described in Sec. 111. In the course 
of this downward evolution, we must split the SM couplings Qi and iu,d,e from their SUSY 
cousins g* and f * ^ ^ as discussed in Paper 1, and also incorporate the distinction between fi 
and ft. As we will see shortly, this complicates the introduction of the electroweak symmetry 
breaking conditions which are incorporated at Msusy- Finally, below the scale mn where 
the two Higgs doublet MSSM transitions to the one Higgs doublet model, we switch from 
the MSSM quark and lepton Yukawa couplings f, to the corresponding SM couplings A,, 
which we reset along with the gauge couplings, after evolving these to Q = Mz {Q = rrit for 
the top quark Yukawa coupling). We then evolve back to the GUT scale where we re-set 
yU = /2 as discussed below, and iterate the solution to the system of RGEs until convergence 
is obtained to the specified precision. 

1. Weak scale boundary conditions 

For the gauge sector, we take as our input at the weak scale the current PDG values [29] 
ttem, «s and sin^ 9w, which are 

a^^{Mz) = 127.925 ±0.016 ; a,(Mz,MS) = 0.1176 ±0.002 ; 
sm'^ew{Mz,ME) = 0.23119 ±0.00014. 

These are the couplings extracted using the effective theory with the electroweak gauge 
bosons as well as the top quark integrated out at Q = Mz- In order to use the SM for 
evolution for Q > Mz, we must match these couplings to those of the full SM, which, to 
two-loop accuracy implies that the SM gauge couplings in the MS scheme are given by 

Although the gauge and Yukawa coupling RGEs along with their tilde eousins form a elosed system, the 
solution to these RGEs may be sensitive to sparticle thresholds through particle decoupling effects. We 
will see that these effects can be especially important in non-universal models where squark mass matrices 
are non-diagonal (in the basis where the corresponding quark Yukawa couplings are diagonal) and the 
squark mass eigenvalues have substantial splitting. 
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Notice that in order to preserve the SU (2) symmetry of the effective theory down to Q = Mz, 
we have, as mentioned earher, integrated out the top quark at Q = Mz rather than at its 
mass as we do for all other particles. This is the origin of the ln(mj/ /i) terms in the matching 
conditions for the gauge couplings above. We emphasize that we decouple all SUSY particles 
as well as the additional Higgs bosons at the scale of their mass: as a result, we do not get 
corresponding jumps in the gauge couplings as these decouple. Our method — which is also 
used in ISAJET — has an important advantage in that it "sums the logs of the ratio of any 
large mass to Mz" , in contrast to the frequently used procedure that uses MSSM evolution 
down to Mz, and then corrects for this via a "single step evolution" (between the heavy 
scale and Mz) to take into account the difference between the running in the MSSM and in 
the SM. 

Next, we convert the values of these gauge couplings in the MS scheme to their corre- 
sponding values in the DR scheme using the relations [32]: 

^ ^ (45a) 

(45b) 

(45c) 



ai(DR) ai(MS) 
1 1 



a2(DR) a2(MS) Qtt 
1 1 1 



a3(DR) a3(MS) 47r' 

and use the results as boundary conditions at Q = Mz when solving the RGEs. 

For the Yukawas, we begin with the quark masses at Q = Mz (the masses of the light 
quarks and leptons at Mz can be found in Ref. [33]), and convert to SM Yukawas using 
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vsM = 248.6/-\/2 as in Ref. [25]. The masses of the first two generations of quarks have 
substantial error, which leads to a corresponding error in their Yukawa couplings. The 
third generation quark masses are more precisely known — we take the top pole mass 
nit = 172 GeV, and rrihlMz) = 2.83 GeV [34]. We rotate the diagonal Yukawa couplings, 
which are in the "quark mass basis" to a current basis using (39) of Paper I at Q = rrit. We 
include SUSY radiative corrections [25] (from IS A JET), with inter-generation quark mixing 
neglected, aX Q = Msusy- Finally, for the numerical results that we present in the rest 
of the paper, we parametrize the Kobayashi-Maskawa matrix as in (50) below, and take 
sin a = 0.2243, sin/5 = 0.0037 and sin 7 = 0.0413, 6/3 = 60° and 6a = = 0. 

2. Electroweak symmetry breaking 

It is traditional to use the observed value of M| to determine the value of /i^ from the 
minimization conditions for the scalar potential in the Higgs sector. The inclusion of thresh- 
old effects causes additional complications for this program. Recall that the electroweak 
symmetry breaking conditions are imposed at Msusy = ^m^^m^^ which is always smaller 
than the mass of the heaviest SUSY particle. As a result, //^ is now conceptually and 
numerically different from which is of course the parameter that enters in the Higgs 
boson potential. Moreover, the Higgs potential depends only on Mfj^ = {rri^j^ + and 
Mfj^ = + so that it is not possible to separate fi^ from the SSB parameters m|^^ 

and Tn^Hi ^^^^ specified at the high scale. Using the combinations, 

(Mi^ + M^J =m^^+m^^ + 2 and 
{Ml, - MlJ = m\ - , 

the tree level minimisation conditions of the Higgs potential can be written as, 

(Mi + MiJ = --1^ (M^ - Mi) - \ {g- + g^) {vl + vj) , (46) 

b = sc (Mi + Mi) . (47) 

The first of these fixes the sum (Mi -|- Mi) in terms of the difference (Mi — Mi). 
Since we know the difference at the GUT scale, we can evolve this down to Msusy (along 
with other SSB parameters) during the iterative process that we use to solve the RGEs. At 
Q = Msusy we use (46) to solve for Mi + Mi, which can be evolved back to the GUT 
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scale. We then use the sum to fix /i = /i at the GUT scale, reset the difference to its input 
value, and iterate. The value of the higgsino parameter /i can then be obtained at all scales 
using (Bl). The 6-parameter usual, be eliminated in favour of tan/3 using (47).^^ 

In our discussion up to this point we have ignored another potential complication 
that arises if m^^ > Msusy- In this case, the heavy particles of the Higgs sector 
decouple, and for Msusy < Q < mn, we only have the light doublet in the effec- 
tive theory that we use to calculate the RGEs. In this case, the heavy Higgs doublet 
mass term [c^ + iP'f) + ■^^ ("^^d + If^f) + sc{b + 6*)] and the Higgs mixing term, 
[sc {m'jj^ + \flf) — sc {rn?u^ + + s^h — c^6*] (and its complex conjugate) in Eq. (35), 
together with tan/3, are frozen at their values aX Q = rriH, while the light doublet mass 
parameter, [s^ + li^f) + ("^1/^ ~'~ l**^!^) " -^c (& + &*)] , along with t^sM = a/w^ + f^, 
continues to evolve to Msusy- The three frozen coefficients together with the evolved mass 
term for the light doublet can now be used to solve for {m'jj^ + \ flf), ("2^^ + \ flf) and the 
complex 6-parameter. We can now find an iterative solution in the same manner as for 
Msusy > mn- 

Before closing this section, we should add that although we have discussed EWSB con- 
ditions only at tree-level, in practice, we minimize the one-loop effective potential including 
effects of third generation Yukawa couplings, but ignoring all flavour-mixing effects in this 
computation. These corrections, which effectively shift the Higgs boson SSB mass squared 
parameters by S„ and S^, respectively, are evaluated by replacing ftM,T in the relations by 
the (3,3) element of the corresponding Yukawa matrices, and with the dimensionful param- 
eters also replaced by the (3,3) element of the corresponding matrix (or the appropriate 

The alert reader will notice that something is amiss in (47): the 6-parametcr is complex, while the right- 
hand-side is manifestly real. The point is that at any one scale, chosen here to be Mgusy, the 6-parameter 
can always be made real. Indeed, the very fact that we have written positive values for the VEVs 
and Vd mandates that b is real and positive at Q = -Msusy- To be specific, we can always make a gauge 
transformation such that just the lower component of the scalar doublet Hu has a VEV, and that this 
VEV is real and positive. Then the minimization of the scalar potential in the Higgs sector requires 
that the VEV of Hd is aligned; i.e. it is also only in its lower component. This alignment is a result 
of the dynamics. Finally, we can redefine the phase of the doublet superfield Hd so that Vd is real and 
positive. This is not compulsory, but is the customary practice that allows us to define tan/3 to be real 
and positive. If instead we write the VEV of the down type Higgs field as \vd\ exp {i9d) and tan/3 = 
the left-hand-side of (47) would have to be amended to bexpiOb- So, although we choose real VEVs and 
a concomitantly real 6-parameter at Q = AfsusYj we retain the complex b in the RGEs, since b will not 
remain real at other scales. 
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frozen value). 

3. SSB parameters at the high scale 

Since our purpose in these papers is to address flavour physics of sparticles in as general 
a way as possible, subject to experimental constraints that seem to suggest that flavour 
physics is largely restricted by the structure of the Yukawa coupling matrices, we thought 
it would be useful to first seek a general parametrization for SSB parameters that does 
not introduce a new source of flavour-violation, but allows for non- universality of model 
parameters. Additional, but uncontrolled, flavour-violation can easily be incorporated by 
allowing for other contributions to the SSB mass and trilinear parameter matrices. We 
use the (s) quark sector to illustrate our arguments, but almost identical considerations will 
apply to (s)leptons, except that in this case we would also have to include additional lepton 
number and lepton-flavour violating matrices in the singlet (s)neutrino sector. 

Within the framework of the -R-parity conserving MSSM, the SSB matrices m^^g and 
^u,d potentially include new sources of flavour-violation, not included in the superpotential 
Yukawa couplings. In order not to introduce a new source of flavour-violation, these SSB 
matrices must be diagonal in the same superfield basis (where the SM fermions and their 
scalar superpartners are rotated by the same matrices) that the superpotential Yukawa 
interactions (renormalized at the same high scale as the SSB parameters) are diagonal. 
Of course, it is impossible to simultaneously diagonalize and f^, but what we mean 
is that the SSB mass matrices that describe the mixing of both left- and right-up type 
squarks, and their trilinear couplings must be diagonal in the basis that the matrix 
is diagonal, and likewise for the down sector. However, since SU{2) symmetry dictates 
that the m? and m^- SSB matrices must be the identical, this doublet squark mass- 
squared matrix must be proportional to the unit matrix, 31, in order to remain diagonal, 
both when the up- or the down-type Yukawa coupling matrix is diagonal. In contrast, 
the matrices mfj and (m|, and a^) can be functions of the up (down) type Yukawa 
coupling matrices (and their Hermitian adjoints) chosen in such a way that these matrices 

This is not equivalent to the requirement that the Yukawa couphig matrix commute with the corresponding 
a-parameter matrix because these non-Hermitean matrices are diagonaUzed by bi-unitary rather than by 
unitary transformations. 
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are simulaneously diagonal when we transform these to the basis where the corrsponding 
superpotential Yukawa couphng matrix is diagonal. 

To find the most general parametrization of the mfj^ and a^ ^ matrices of the type that 
we are looking for, we first note that these, respectively, transform in the same way as the 
matrices (fj^fj^rf)" and fu^d (^fldfw.dj ! where n is any integer. Thus any linear combination 
of these matrices (with ?T, = 0,l,2---)is guaranteed to be diagonal in the basis that iu.d is 
diagonal (at the high scale at which we input the SSB parameters). The only question, then, 
is just how many terms we need to allow in the linear combinations that make up or 
SLu^d, to guarantee the most general form for these SSB matrices, so that fiavour violation 
enters only through the superpotential Yukawa coupling matrices? This is easiest to see in 
the diagonal basis for the Yukawa couplings. The SSB matrices are also diagonal in this 
basis, and so are completely specified by Ug diagonal elements, where Ug is the number of 
generations. Transforming to a general basis does not alter the number of parameters that 
we need: we thus know that we must have Ug terms in each of the linear combinations for 
m^^ and for SLu^d that we discussed above. For the MSSM with Ug = 3 generations, we thus 
parametrize the SSB sfermion mass and a-parameter matrices at the high scale as, 

^l,L = mlQ^L}o^ + TQ^L, (48a) 

™V,D,E = ''^Iu,D,E}o['^U,D,e'^ + Ru,D,E^I,d,e^u,d,e + Su,D,EifI,d,eK,d,e)'^] + T^U,D,E , (48b) 
d,e[^{w,d,e}olJ- + Wu,d,e^u4,e^u,d,e + Xu4,e{S-^^d,e^u,d,e) ] + 'Ziu4,e , (48c) 

where iu,d,e are the superpotential Yukawa coupling matrices in an arbitrary current basis at 
the same scale at which the SSB parameters of the model are specified. Here, cu,d,e = or 
1 is introduced only to allow the facility to "switch off" the universal term if desired. The 
matrices 'Tq,l,u,d,e and Tiu^^e have been introduced only to allow for additional sources of 
fiavour-violation not contained in the Yukawa couplings. Setting Tq 2,,(7,d,£; = '^u,d,e = 
gives us the most generation parametrization of the three-generation i?-parity conserving 
MSSM where the superpotential Yukawa interactions are the sole source of fiavour viola- 
tion. In other words, any top-down theory of fiavour, will specify the form of the matrices 
Tq,l,u,d,e and Zu,d,e along with the coefficients m^Q, A,o, c,, R,, S,, W, and X, that appear 
above. 

The physics behind our anstaz (48) for SSB parameters, with T, = Z, = 0, is that we 
assume that the SUSY breaking mechanism, for reasons that are not understood today, does 
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not introduce a new source of flavour violation.^ We stress that this differs from the minimal 
flavour violation [36] assumption, where the new physics may introduce an additional flavour 
violation in a controlled way. Our ansatz is thus a special case of the minimal flavour 
violation scenario. Of course, the form of (48) is not invariant under renormalization group 
evolution, and at the weak scale we will end up with the same form for the SSB parameters 
as the minimal flavour violation ansatz (suitably generalized from the flrst paper of Ref. [36] 
to include higher order terms in the Yukawa couplings), but with appropriate relations 
between the coefficients of the various terms in the minimal flavour violation formulae for 
SSB parameters. We will return to the difference between minimal flavour violation and 
(48) in Sec. VI. 

We have created a code (to be incorporated into ISAJET) that allows the reader to 
choose the arbitrary values for the various coefficients as well as the matrices Tq^l,u,d,e and 
Zu,d,e that appear in (48) as boundary conditions for the RGEs for the SSB parameters of 
the MSSM. This code provides a tool for a study of (s)quark flavour violation in sparticle 
processes in a arbitrary theory that reduces to the MSSM at the scale where the SSB 
parameters are specifled. We can evolve the gauge and Yukawa couplings from the weak 
scale to this high scale, use the boundary conditions (48) as inputs for the SSB parameters, 
and proceed to determine the iterative solution to the RGEs as described above. The 
familiar universal mSUGRA boundary conditions are reproduced by setting cj/^.s = 1; 

^{Q,L}0 ~ ''^{U,D,E}0 ~ ^'o'l ^{M,d,e}0 = ^O'l Ru,D,E = Su,D,E = Wu4,e = Xu,d,e = 0; Tq^l = 

Tu,D,E = Zu,d,e = in (48). 

A. Quark Yukawa Couplings 

We begin the discussion of our numerical results by showing the magnitude of the complex 
elements of the up quark Yukawa coupling matrix above Q = mn, and A„/ sin/5 below 
Q = niH, for the mSUGRA model with mo = 200 GeV, mi/2 = -400 GeV, Aq = -200 GeV, 
tan/5 = 10 and /i > in Fig. 1.^^ The approximate spectrum for this illustrative scenario is 

For a complementary approach, where supcrsymmetry breaking is the origin of fermion flavour, see 
Ref. [35] , and the extensive bibliography therein. 

The reader may well wonder why we choose the gaugino mass to be negative and of the opposite sign 
to that of the /i parameter, a relative sign that is "apparently disfavoured" by the Brookhavcn — 2 
measurement. The reason is that our convention for the sign of fi [4] (as well as the one used in ISAJET) 
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mSUGRA: m =200 GeV, m =-400 GeV, A =-200 GeV, tanP=10, |X>0 
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FIG. 1: Evolution of the magnitudes of the complex elements of the up-quark Yukawa coupling 
matrix of the mSUGRA model with mo = 200 GeV, mi/2 = -400 GeV, Aq = -200 GeV, tan/3 = 
10 and /i > in the basis where this matrix is diagonal at Q = m^. For Q > tjih (— 631 GeV) 
we plot |(fu)ij| whereas for Q < rriH, where the effective theory includes just one scalar Higgs 
doublet, we plot KAm)^! / sin/3 which is equal to |(fu)ij| at Q = mu- In all the figures we take 
mt = 172 GeV. 

shown in Table I. Fig. 1 should be compared with Fig. 1 of Paper I, where we had shown 
the evolution of the same Yukawa coupling, but in a simplified scenario where the SUSY 
thresholds were clustered in two regions, one at 600 GeV and the other at 2 TeV. 

The most striking feature of the figure are the dips in the off-diagonal elements of the 

is opposite the one that is usually used. To compensate for this, ISAJET internally (and unknown to the 
user) flips the sign of the gaugino masses relative to what the user uses as an input. Since the physics 
depends only on the relative sign between the gaugino masses and ^, this sign flip is equivalent to a 
flipped sign for fi (together with sign flips for Aq and 6), making it appear that the ISAJET output is 
in accord with the usual convention for the sign of ^. We should also mention that here and in Ref. [4], 
the convention for the sign of the a is opposite that to the one in Ref. [18] and of ISAJET. To reproduce 
the results with the sign conventions for mSUGRA parameters that we use in this paper (and in Ref. [4]) 
using ISAJET, we should merely switch the sign of the gaugino masses in the ISAJET inputs. Our results 
(except for inter-generation mixing effects) can thus be obtained by running ISAJET with mo ~ 200 GeV, 
mi/2 = 400 GeV, Aq = -200 GeV tan/? = 10 and > 0. 
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TABLE I: The location of the thresholds for our canonical mSUGRA case in Fig. 1 and in several 
subsequent figures. 

iu matrix. The aligned dips at Q ~ 650 GeV, common to all the off-diagonal elements, 
occur because of the change in the sign of the coefficient of the f^fj-type terms that drive 
the growth of these off-diagonal elements from zero at Q = rrit. This has been discussed in 
detail in Paper I, where we have also explained why the magnitudes of all the off-diagonal 
Yukawa couplings vanish at essentially a common value of Q. The presence of the second zero 
at the higher value of Q, just in the (2,3) and (1,3) elements, is accidental. It occurs because 
of conspiracies between terms in the corresponding /3-function as the left-type squarks are 
decoupled. Notice that the lowest four curves, though they do not have this additional dip, 
show kinks at these same values of Q, corresponding to the decoupling of these squarks. For 
several other mSUGRA cases, we have checked that while squark decoupling causes kinks 
in the curves, the coupling does not drop to even close to zero for a second time, in contrast 
to the behaviour in our illustrative example in the figure. 

The evolution of the down-quark Yukawa coupling matrix for the same mSUGRA point is 
shown in Fig. 2. We have shown these couplings both in what we will subsequently refer to 
as our "standard" current basis where the up-type Yukawa couplings are diagonal at Q = rrit 
in the left frame, and in the basis where the down-type Yukawa coupling matrix is diagonal 
at Q = rrit in the right frame. The matrices in the two bases are connected by the KM 
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mSUGRA: m =200 GeV, m =-400 GeV, A =-200 GeV, tan|3=10, |^>0 
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(1,3) _ 
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(1.2) _ 

(2,1) - 



10000 le+08 le+12 le+16 10000 le+08 le+12 le+16 

Q (GeV) Q (GeV) 

FIG. 2: Evolution of the magnitudes of the complex elements of the down-quark Yukawa coupling 
matrix of the mSUGRA model with tuq = 200 GeV, ?ni/2 = -400 GeV, = -200 GeV, tan/3 = 
10 and /i > 0. In the left-frame, we show the elements of the matrix in the same basis as in 
Fig. 1, where A„ is diagonal at Q = nit, whereas in the right frame, we show these elements in the 
basis where the is diagonal at Q = mt. For Q > mu (— 631 GeV) we plot |(f(i)ij| whereas for 
Q < rriH, where the effective theory includes just one scalar Higgs doublet, we plot |(A(i)jj| / cos/3 
which is equal to | (frf)jj- 1 at Q = mn- 

matrix, in an analogous way to (39). The curves in the left frame are all smooth (except for 
the small kink in the curves for | (fd)j3 | that occurs because of the SUSY correction to the 
bottom quark Yukawa coupling [25]), and do not show the dip to zero that appeared in the 
previous figure. This is not surprising because underlying the explanation of this dip was 
the fact that the off-diagonal elements evolved from zero at Q = [14]. Notice also that in 
this frame the off-diagonal elements are not necessarily smaller than the diagonal elements 
even for Q ~ 1 TeV. The magnitudes of the off-diagonal matrix elements in the frame on 
the right, which do start at zero aX Q = rrit show the anticipated aligned dips, except that 
the location of the dip is shifted considerably to the right, relative to Fig. 1. This shift is 
not difficult to understand. The large top quark Yukawa coupling in the SM governs the 
evolution of the off-diagonal elements of the down-type Yukawa couplings, causing them to 
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evolve much more rapidly from zero in the right-hand frame of Fig. 2, so that in order to 
evolve back to zero after the sign flip in the /^-function due to the additional Higgs and SUSY 
particles, a longer evolution distance is needed in the present case. Furthermore, beyond 
the Higgs boson threshold the off-diagonal elements of in Fig. 1 are accelerated to zero 
on account of the fact that the down-type Yukawa couplings f^, that enter in the evolution 
of these elements, are enhanced by a factor ~ 1/ cos/3, pushing the dip in this figure to a 
low value of Q. 

We now turn to briefly discuss what happens when we allow non-universality of GUT 
scale squark mass parameters, so as to split the squarks more than in mSUGRA. Recall 
that if the squarks all decouple together, all that happens is a change in the slope of the 
/3-function. If, however, the squark masses are not all the same, our decoupling procedure 
entails an additional rotation to the squark mass basis. If this basis differs significantly from 
our "standard" basis in which the quark Yukawa coupling matrix is diagonal (at Q = rrit), 
one may expect considerable deviation in the evolution of the off-diagonal elements from 
the mSUGRA case for Q values in between the highest and lowest squark thresholds. The 
introduction of non-universal squark mass parameters via non-vanishing values of R, and 
S, in (48) never leads to significant effects because the rotation from the "standard" basis 
to the squark mass basis is small by construction. 

The question then is whether we can have large deviations from Fig. 1 and 2 via the 
T, (or Z.) matrices. To examine this, we set all GUT scale inputs to be the same as the 
mSUGRA case in Fig. 1 except that we now take, 

m'^u,D}o = 5 (49a) 
Tu,D = diag {10000, 40000, 90000} GeV^ . (49b) 

We need, of course, to specify the basis in which the squark mass matrix is diagonal. If 
this is the "standard" current basis, we have checked that although the right squark masses 
are now significantly split relative to the mSUGRA case, except for detailed changes in the 
evolution for Q values in between the squark thresholds {e.g. the dip at the higher value of 
Q in the upper curves in Fig. 1 develops more structure, and there are jumps in the (1,2) 
and (2,1) elements), the evolution is not altered in any important way at large values of 
Q. This is because the rotation between the "standard" basis and the squark mass basis is 
again small. 
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If we take instead the right-squark matrices in (49) to be diagonal in a completely different 
basis specified by unitary matrices Vl(u), y^R^u) and y^R^d) that take us to our "standard" 
basis, and specified to be of the form, 



V 



(50) 



(where Sa = sin a, Cq, = cos a, etc.) with randomly chosen values of a, P, 7 and 5/3 (we take 
= ^7 = for simplicity), we may expect the evolution of Yukawa couplings to depart 
from the corresponding evolution in mSUGRA.^^ We emphasize that this scenario (which is 
only of pedagogical interest) where are diagonal in the basis where the corresponding 

Yukawa coupling matrices have large off-diagonal elements includes potentially very large 
flavour-violation in the singlet squark SSB mass matrices, and is likely excluded. 

In Fig. 3 we illustrate the evolution of the magnitudes of the elements of the up-quark 
Yukawa coupling matrix with, 

Vl{u) : a = 2.053, p = 0.254, 7 = 2.030, S/3 = 0.4829 , (51a) 
Vr{u) : a = 1.188, p = 2.218, 7 = 0.763, S/3 = 0.87 , (51b) 
VR{d) : a = 1.904, p = 2.947, 7 = 1.847, 5/3 = 1.14 . (51c) 

Note that, just as in Fig. 1, we have plotted these elements in our "standard" basis where the 
Yukawa coupling matrix is diagonal at Q = mj. The GUT scale matrices and mj^ will 
be random-looking Hermitian matrices in our "standard" basis. We see that the evolution of 
the diagonal elements is not significantly altered from mSUGRA. This is because although 
there is a large mismatch between the squark mass basis and our "standard" basis, this 
mismatch is operative over the small range of Q between the highest and lowest thresholds 
and so has little impact on the largest elements. The two largest off-diagonal entries (i.e. 
(fu)23 and (fu)i3) similarly do not change significantly from Fig. 1, but all other entries are 



The matrices Vi^ii(it) should be numericaUy unitary to high accuracy. Otherwise, numerical errors from 
inverting the up Yukawa couplings from this new current basis to our "standard" basis will leave residual 
off-diagonal elements rather than zero even at Q = mt . If the size of these elements is comparable to the 
values of the smallest off-diagonal elements at values of Q substantially away from mt, it is clear that 
our solutions will be dominated by the error from the non-unitarity of the 'VL,Riu) matrices. A similar 
consideration apphes to Vij(d). 
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mg=200 GeV, mj/2=-400 GeV, Aq=-200 GeV, tanP=10, |i>0 

™%Dio=0' TuD= diag( 10000,40000,90000) GeV^ and Vj^(u,d)#l 
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FIG. 3: The same as Fig. 1 except that the GUT scale right-squark mass matrices are now given 
by (49) in the basis specified by (51) of the text. Note, however, that these matrix elements are 
shown in the "standard" current basis where the matrix is diagonal at Q = mt. 

greatly altered. While it may seem that the values of these matrix elements at a large scale 
is quite irrelevant phenomenologically, the altered form of the Yukawa coupling matrix at 
the high scale could be of relevance to model-builders. 

We have seen that the magnitudes of the off-diagonal elements in Fig. 3 remain small 
because the splitting in the squark spectrum is limited to (9(100 — 1000) GeV. The natural 
question then is whether we can get these to be larger by choosing extreme intra-generational 
squark splitting. Even putting aside potentially unacceptable flavour- changing effects that 
might result, this is not easy. In general, such a GUT scale splitting also has a large value 
for S = rnjj^ —m'jj^ + Tr [nig — m| — 2m^ + + m|;] , which pulls the other squarks also 
to large masses, so the squark mass splitting is reduced by RGE effects. It may be possible 
to obtain split squarks by adjusting S to be zero {S is then invariant under renormalization 
group evolution), but we have not investigated this here. 
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B. Gaugino Mass Parameters in Split SUSY Models 

We have seen in Sec. IV A that the evolution of the electroweak gaugino mass parame- 
ters acquires a dependence on the /i-parameter and vice-versa, once threshold effects from 
splitting in the Higgs boson sector are included. In models where rriH ^ |Afi,2|; I^i2l> 
the effect of the term explicitly dependent on /i in (B2)-(B5) (and the corresponding terms 
dependent on the gaugino mass parameters in the RGE for fi) may be significant, so that 
the relation M2/M1 = a^joix expected in many models is modified. We should keep in mind 
that two loop terms will, in general, also alter this relation. Our point is that we should 
expect threshold corrections, from the [i term in the RGE, as well as from the decoupling 
of sfermions, to be comparable to (or even larger than) the two loop modifications, and so 
need to be included in a quantitative analysis. Within the mSUGRA context, we have small 
values of and hence, mn ^ in the hyperbolic branch/focus point (HB/FP) region 
[37] which occurs for large values of mo, and is one of the regions selected out [38] by the 
relic-density measurement [39]. We mention here that the location of this HB/FP region is 
significantly altered by the inclusion of the threshold corrections. 

These considerations led us to examine the evolution of the gauge couplings and the 
electroweak gaugino parameters for a relic-density-consistent mSUGRA model point in the 
HB/FP region with mo = 3075 GeV, mi/2 = -600 GeV, Aq = 0, tan/5 = 10 and /i > 0, 
for which {/j, M^, M2) ^ (306,-258,-496) GeV at the weak scale.^'^ We found that at 
the two loop level with all threshold effects included, M2/M1 = 1.919 to be compared to 
M2/M1 = 1.881 obtained without including threshold effects. Thus, although threshold 
effects actually bring us closer to 02/0^1 = 1.964, their inclusion is clearly necessary for a 
quantitative analysis of mass parameters that may be extracted at an e~^e~ linear collider, 
where a precision of better than 1% will be possible if charginos are kinematically accessible. 

The threshold corrections to gaugino mass parameters can be much larger in the so-called 
split SUSY model [40] that has received considerable attention in the recent literature. In 
this scenario, the naturalness of the scalar Higgs sector (which we view as one of the primary 
motivations for weak scale SUSY) is abandoned, while gauge coupling unification and the 
neutralino dark matter candidate of i?-parity violating models are preserved. Gaugino mass 

Without threshold corrections, ISAJET gives a similar spectrum for mo ~ 3660 GeV. 
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Split Supersymmetry 
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FIG. 4: The evolution of the gaugino mass ratio M2/M1 at the one (dashed) and two (solid) loop 
levels, along with the two-loop evolution of a2/ai, for a split SUSY model where scalar masses 
are around 10^ GeV, with gauginos and higgsinos at the weak scale. Gaugino mass unification 
is assumed and as everywhere else in this paper, the gaugino mass parameters are negative. The 
other parameters are as mentioned in the text. 

parameters and are assumed to be at the weak scale, while scalar mass parameters are 
set to be at an intermediate scale. This means that sfermion masses as well as niH are very 
large (with the SM Higgs doublet fine-tuned to be light), so that charginos and neutralinos 
are the only new particles (other than a SM Higgs boson) at the weak scale. 

As an illustration, in Fig. 4 we plot the variation of the ratio M2/M1 at the one- loop 
level (dashed), as well as at the two-loop level (solid), with the renormalization scale Q, 
along with the two-loop value of a2/cti- We show results, first where the value of is 
set exactly in between \Mi\ and IM2I so that the lightest neutralino acquires a significant 
higgsino component, to qualitatively mimic mixed higgsino dark matter. We have checked 

In the absence of a real theory of split SUSY, we should view this figure only as a qualitative illustration 
of potentially large threshold effects. Here, we take the sfermion mass parameters to be 10^ GeV at 
Q = Mgut, = —350 GeV and ^0 = 0. Since it is not possible to satisfy the EWSB conditions 

except when tan/3 is hierarchically large — this would cause down-type Yukawa couplings to become 
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that the M2/M1 values do not change in a significant way for yet larger values of tan/?. It 
is clear that the relation ^ = ^ is violated at the several percent level by the threshold 
corrections, without which the M2/M1 lines would have continued with the same slope that 
they have above 10^ GeV to low values. For this point, the 2sc x /i term that explicitly 
appears in the RGE is very small so that the result is independent of the sign of /i: most 
of the difference is an effect of the sfermion loop contributions being switched off below 
10^ GeV. To gain some idea of how large the effect of this /i term might be, we have also 
shown M2/M1 for /i = —2 TeV, with sign{fiM2) > and tan/3 = 4. Since such a large value 
of /i would be totally incompatible with the measured relic density and small values of tan (3 
are unnatural in these models without some modification to the EWSB sector, the reader 
should view these curves only as a guide to how much the gaugino mass ratio may deviate 
from its "unification value" . The difference between the two cases is almost entirely due to 
the different choice of /i. We see that the gaugino mass unification condition will, in this case, 
be violated by ~ 10%. If instead we choose the opposite sign for M2fi, but keep = 2 TeV, 
this "large \fi\ line" would be lower than the corresponding line with fi = (Mi + M2)/2 by 
about the same amount that it is higher than this line in the figure. Clearly, increasing the 
splitting between the scalar and the gaugino/higgsino sector of the theory will cause even 
further violation of the unification condition, and ~ 20% effects appear to be plausible if 
the scalars are instead at the 10^^ GeV scale. 

C. Trilinear Couplings 

We begin our discussion of the trilinear couplings by returning to the mSUGRA case 
considered in Figs. 1 and 2. The magnitudes of the individual entries of the trilinear coupling 
matrices a„ and are shown in Fig. 5, where we plot (a) | {slu)^^ \ and (b) \ (ad)-^- | in our 
"standard" current basis where the up-type quark Yukawa coupling matrix is diagonal at 
Q = nit, along with | {3id)y \ in the basis where the down-type quark Yukawa coupling matrix 
is diagonal a.t Q = rrit in frame (c). We terminate the curves a.t Q = rriH since below this 

non-pcrturbative — we treat /i and tan/3 as phenomenological parameters, and fix mn to be 10^ GeV 
in this figure. The parameters mj^^ and mj^^ (indeed all scalar mass parameters) are never needed since 
the RGEs for gaugino masses, ^ and the a-parameters, and the dimensionless couplings, form a closed set 
even at the two-loop level. Sfermion masses only enter via the location of thresholds. 
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FIG. 5: The magnitude of the elements of the trihnear couphng matrix a, for the mSUGRA model 
in Fig. 1. We show (a) | {siu)ij |, (b) \ {sid)ij !> in the basis where up-quark Yukawa couplings are 
diagonal at Q = mt, and (c) | {sid)ij \ in the basis where down-quark Yukawa couplings are diagonal 
at Q = nit- The curves extend between Q = mn and Q = Mqut- 

scale we have a single Higgs scalar doublet model, and the trilinear couplings evolve only 
as part of a linear combination with /i*f^^''' as discussed in Sec. IV B. We see that the curves 
in frame (a) show a simple dip structure indicating that the real and imaginary parts of a^^^ 
are really monotonic functions of Q that pass through zero together, in a manner similar to 
the elements of the Yukawa coupling matrix [14]. The actual location of the zero is somewhat 
harder to analyse because even though a„ obtains off-diagonal components only because the 
down-quark Yukawa matrix is not diagonal a.t Q = rrit, the matrix a^ is off-diagonal even 
at the GUT scale. The off-diagonal elements of a^ in frame (b) start off with a much larger 
magnitude in the "standard" Yukawa basis at Q = Mqut because the corresponding Yukawa 
coupling matrix has large off-diagonal pieces. In this case, the evolution of these off-diagonal 
elements receives significant contributions from all entries in the RGE (unlike the evolution 
of the off-diagonal elements in frame (a) or of the off-diagonal Yukawa couplings discussed 
at length in Paper I where contributions from the off-diagonal down-type Yukawa matrices 
govern the evolution), and never go through zero; the situation is similar to that in the first 
frame of Fig. 2. We see from frame (c) that the magnitudes of (a^^)^^., in the basis that the 
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down-type Yukawa coupling matrix is diagonal at Q = rrit, again show the characteristic dip 
structure indicating that the off-diagonal elements increase in magnitude from their value 
aX Q = Mqut to some maximum magnitude at an intermediate scale, but then smoothly 
reverse direction and thereafter evolve monotonically through zero to the low scale. The 
elements in frames (6) and (c) are, of course related by, (39). We note that because the 
off-diagonal elements in frame (c), are driven by the larger (and significantly off-diagonal) 
"up- type" Yukawa couplings, these are bigger than those in frame (a) where it is the down- 
type Yukawa coupling matrix that has the significant off-diagonal elements, and so largely 
determines the entries. 

Finally, in Fig. 6 we consider a model with non-universal values of a-parameters, but 
where these are not a new source of flavour violation. Specifically, we consider a model with 
the same values for mSUGRA parameters as in Fig. 5, but with non-zero values for W and 
X in (48c) (with A{„rf}o = ^o); and illustrate the evolution of (a) Re(a„)23 and (6) Re(a„)32. 
We have checked that the imaginary parts of these matrix elements are about four orders of 
magnitude smaller. The striking feature of frame (a) is that the various curves which start 
with very different values of Re(a„)23 at Q = Mqut, appear to focus to a common value 
at the low scale. We have checked, however, that although they all cross at Q ^ 1.5 TeV, 
they do not all converge at precisely the same value of Q. This apparent convergence, which 
persists for other values of mSUGRA parameters, is sensitively dependent on the special 
GUT scale boundary conditions for a^ that we have used. We have checked that if instead 
we use a general matrix Zu in (48c), the corresponding evolution is completely different. We 
do not have a good explanation for the seeming convergence in frame (a), and only note 
that it is not generic to all elements of a^ as evidenced, for example, by the corresponding 
evolution of Re(au)32 in frame (b) of the figure. 

D. Soft Masses 

We begin our discussion of the evolution of the scalar mass SSB parameters by showing 
in Fig. 7 the evolution of the magnitudes of mfj in our "standard" current basis for the 
mSUGRA model with the same parameters as in Fig. 1. The diagonal matrix elements 
start from a common value mo and increase as we go to the weak scale because of gauge 
(and gaugino) interactions. The splitting between the (1,1) and (2,2) elements that occurs 
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mQ=200 GeV, mj/2=-400 GeV, Aq=-200 GeV, tan|3=10, |i>0 
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FIG. 6: The evolution of (a) Re(au)23 and (6) Yie{a^)^2 for a model with the GUT scale values of 
a-parameters set as in (48c) for several values of W and X shown in the legend (with A{u^d}o = ^o)- 
The imaginary parts of the matrix elements are about four orders of magnitide smaller. The curves 
are in the same order as in the legend. The GUT scale SSB scalar and gaugino mass parameters 
are assumed to be universal, with the same value as in Fig. 5. We terminate the curves at Q = mn 
at the lower end. 

because of the Yukawa couplings is too small to be visible in the figure. The (3,3) element is, 
however, reduced significantly on account of the large (3,3) entry in f^- Notice that the curves 
become flat once the squarks are all decoupled. The magnitudes of the three independent 
off-diagonal elements of the Hermitian matrix mfj start from zero at Q = Mgut, and rapidly 
rise because has off-diagonal entries at Q = Mqut^ the much more off-diagonal matrix 
affects the evolution of only at the two-loop level. Note the break in the vertical scale 
in the figure. The ordering of the magnitudes of the off-diagonal elements of mfj can be 
simply gauged from the up-type Yukawa coupling matrix. These off-diagonal elements start 
from zero at Q = Mqut, evolve to a maximum magnitude, then smoothly reverse direction 
at an intermediate scale and then continue to evolve monotonically all the way to the weak 
scale. The dips in the figure occur where the real, and simultaneously the imaginary, part 
of (m^) ■ ■ changes sign during the course of its evolution to Q = rrit where we terminate the 
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mSUGRA: m =200 GeV, m =-400 GeV, A =-200 GeV, tanP=10, ^>0 
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FIG. 7: The scale dependence of the magnitudes of the entries of the matrix for the mSUGRA 
model with parameters as in Fig. 1 in the basis where the up-type quark Yukawa coupling matrix 
is diagonal at Q = m^. Note the we have broken the vertical scale at 0.1 GeV^ and also used a 
different scale above this to better show the splitting of the (3,3) entry from the other diagonal 
entries. 

plot. 

Before proceeding further we draw the reader's attention to a technical point that is 
important for the numerical solution of the RGEs, once squark decoupling is included as 
described in Sec. Ill B. As we have already explained, in order to correctly implement the 
decoupling procedure for values of Q below the highest sfermion threshold, at each step 
we need to rotate to the basis where the SSB squark mass squared matrices are diagonal, 
which, in turn, requires us to obtain the unitary matrix R that relates our "standard" basis 
to this "squark mass basis". After evaluating the right-hand side of the RGEs, we always 
rotate back to our "standard" basis using "R) . While this is straightforward in principle, the 
practical problem is that when two squark eigenvalues become very close — this is always the 
case in mSUGRA because the up and charm squark mass parameters only evolve differently 
because of effects of the small first and second generation quark Yukawa couplings — RR^ 
develops non-zero off-diagonal matrix elements due to numerical noise at the 10" ^° level (The 
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mSUGRA: mg=200 GeV, mj^2=-400 GeV, Ap=-200 GeV, tanP=10, n>0 

"standard" current basis '^^^^''^ ^^ere fjm) diagonal 
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FIG. 8: The scale dependence of the magnitudes of the entries of the matrix uIq for the mSUGRA 
model with parameters as in Fig. 1 in the basis where the up-type (left frame), or the down-type 
(right frame), quark Yukawa coupling matrix is diagonal aX Q = mt- Note the we have broken the 
vertical scales to better display the matrix elements. 

noise level depends on the computer system.). This then ruins the delicate cancellations that 
are necessary to obtain the tiny magnitude of (m^)-^2 seen in the bottom curve of Fig. 7. 
Our procedure for dealing with this is detailed in Appendix A. It is for essentially the same 
reason that we had to use the manifestly unitary form (as in (50)) for the matrices Vi('u), 
V^(u) and \ ji{d) when we discussed the evolution of Yukawa couplings when the matrices 
m^£, were diagonal in a general current basis in Fig. 3; see the footnote just after Eq. (50). 

Fig. 8 shows the evolution of itiq for the same mSUGRA point. In the left frame, we 
show the magnitudes of the elements in our "standard" basis, where the up quark Yukawa 
coupling matrix is diagonal at Q = fnt. The frame on the right shows the magnitudes 
of the elements of this same matrix, but in the basis where the down-type quark Yukawa 
coupling matrix is diagonal at m^. We have checked that the large difference in the size 
of the off-diagonal elements in the two frames is indeed accounted for by the fact that the 
corresponding matrices are related by (40). Unlike in Fig. 7, there is no dip in the magnitudes 
of the off-diagonal elements because they evolve monotonically from zero at the GUT scale. 
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FIG. 9: The scale dependence of (111^)23 for several sets of values of Ru,D and Su,d that appear 
in (48b) for the model with non- universal GUT scale SSB squark mass parameters, but cu,d = 1- 
The curves are in the same order as the legend, and all the other parameters are set as in Fig. 7. 

until the squarks are all decoupled. 

To understand how the non-universal boundary conditions in (48b) impact the evolution 
of the squark mass matrices we examine a model with non- zero values of Ru,d and Su,d, 
but with universal gaugino masses and a-parameters. As an illustration we show in Fig. 9 
the value of (1x1^)23 for the set of values of Ru,d and Su,d shown in the figure, with all other 
parameters set as in Fig. 7 (including cu^d = 1) so that Ru,d = Su,d = corresponds to 
mSUGRA model in this figure. We see that with non-zero values of Ru,d and Su,d, {''^u)23 
already starts off with a substantial value (positive or negative) aX Q = Mqut, and evolves 
slowly with Q. This situation is qualitatively similar to that in Fig. 7, once (m^)^^- has 
evolved away from its value at Mqut? and has had a chance to grow from zero. However, 
because the curves start of with rather large values at the GUT scale (except for the middle 
mSUGRA curve) the evolution does not take them through zero for any value of Q > Trit. 
As a result, the dip which was the most prominent feature of Fig. 7 is absent, except in the 
middle curve which does cross zero for Q ~ 2.5 x 10^ GeV. 

Finally, in Fig. 10 we return to the non-mSUGRA case that we considered in Fig. 3. 
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FIG. 10: (a) The same as Fig. 7 except for the non-universal model considered in Fig. 3, and (b) the 
magnitude of the corresponding elements of m|, for the same scenario. The elements are plotted in 
our "standard" current basis but are exactly the same in the basis where the down quark Yukawa 
couplings are diagonal. 

We emphasize again that in this case mfj and m|, are diagonal only in the basis where the 
superpotential Yukawa coupling matrices have large off-diagonal elements so that we would 
expect that this model includes new and potentially large sources of flavour-violation in the 
singlet squark SSB mass matrices that may well be excluded by data. As in Fig. 7, we show 



mfj)^j I in the basis where the up-quark Yukawa coupling matrix is diagonal at Q 



rrit. 



We see from the figure that the matrix mfj has large off-diagonal entries at Q = Mqut- As 
expected, the gauge (and gaugino) interactions cause the diagonal entries to rapidly increase, 
whereas the off-diagonal elements which do not "feel these terms" evolve much more slowly 
with Q. Note also that because (m^)^-,^ < (111^)22 at (5 = Mqut, Yukawa coupling effects 
draw them closer as we go to low scales. We also remark that the negative GUT scale value 
of S tends to reduce the diagonal elements of as we go to low scales but pulls up the 
corresponding elements of m|,. Although the renormalization group evolution increases the 
gap between the off-diagonal and diagonal elements at the low scale, notice that the off- 
diagonal elements are separated by just one order of magnitude from the difference between 
the diagonal elements, so that we may expect large flavour mixing between the SU (2) singlet 
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up squarks. This mixing, if anything, is even larger in the down squark sector as can be seen 
in frame (b). A careful evaluation of inter-generation squark mixing is clearly necessary for 
any discussion of flavour-violation in squark decays, the subject of the next section. 



VI. FLAVOUR-CHANGING SQUARK DECAY: AN ILLUSTRATIVE APPLICA- 
TION 

We are now ready to apply the considerations of our analysis to flavour violating de- 
cays of squarks. In models with additional sources of flavour violation in the SSB sector 
through non- vanishing values of the matrices T, and/or Z, in (48), at least some the SSB 
mass and a-parameter matrices will not be diagonal in the basis that quarks are diagonal, 
and we will have flavour- violating decays of squarks at the "tree-level", i.e. even without 
any renormalization group evolution of the SSB parameters. If the elements of these ma- 
trices take arbitrary values comparable to MgusY? flavour- violating quark-squark-neutralino 
and quark-squark-gluino couplings become comparable in magnitude to the corresponding 
flavour-conserving vertices, resulting in the well-known SUSY flavour problem. Minimal 
flavour violation [36] has been suggested as a way of introducing flavour violation into the 
SSB parameters in a controlled way. We present formulae for flavour violating couplings 
and partial widths for the corresponding decays of squarks, and then go on to a numerical 
analysis of these decays. For definiteness, we will focus on up-type squarks. 

A. Squark decay width 

In the current basis the up-type squark mass terms in the Lagrangian are. 




(52) 



where the (6 x 6) up-squark mass matrix takes the form. 




(53) 



Here, the indices l,m label left-handed squarks {uu = [uL,CL,tL)), while r,s label right- 
handed squarks {uRr = (m_r, c/{, t/j)). The elements of the squark mass matrix are given 
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by 



Im 



l{i'Jl),m+[^-'-i){<-^''i)^^rn, (54a) 

^fJf:)..-y K-^^)'^.., (54b) 

(ML),^ = -v^ (aj;, + {f,*it)l . (54c) 

Of course, if we are evaluating (A4|^)^^ for Q < mu, it is the combination 
[—s (a„)^^ + c (Ai*^^")^^] (for which the RGE is (B21) in Appendix B) that enters the mass 
matrix. 

The physical squarks {u^^ , cr = 1 — 6) are obtained by diagonalizing the squark matrix 
(53) by a unitary transformation lA, and are given in terms of the current squarks by, 



[^U)rs + 




(55) 



where, for later convenience, we write the (6 x 6) matrix lA in terms of two (3 x 6) blocks Ul 
and IA.R. We label the physical squarks in order of their mass, with being the lightest, 
and being the heaviest. Similar considerations hold for down- type squarks. 

The couplings of the physical squarks to neutralinos can readily be evaluated by inserting 
the squark mass eigenstate fields obtained by inverting (55) into the obvious generalization. 



(qRsqaZi) =qR,Zi 



Pr, 



qa + h.c. 
qa + h-c. 



(56a) 
(56b) 



q Jsa^q 

of the quark-squark-neutralino couplings in Ref. [4] that incorporates the fact that the matrix 
couplings g and f are generally different from the usual gauge and Yukawa couplings. Here, 
we have used (f^^)sa, ^q)ma and v'f' to signify (C«)sa, (f^)ma and vf for q = u. We have 
the same form for the interactions of down quarks and squarks, but we must then remember 



to use Vq^ = . The matrices and that appear in the couplings generalize to. 



B 



V2 



V2 



3^2 



B 



z. 



3V2 



9.-1 ,« 
^4 ) 



ma (i) 



'ma (i) 
^4 



(57a) 

(57b) 
(57c) 
(57d) 
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in direct correspondence with (8.87) of Ref. [4]. The Lagrangian for physical squarks then 
takes the form, 



C {q^qaZ?j 3 t'^Z, Pl + (/3|J^^ Pr] Qa + h.c. 



(58) 



with 



oc 



/3" 



Ur)L. ( Al 



am \ Zi 



/3 



z^ 

am K'-^Zi 



E2(X>i)L { A' 



(-OS;«(W^)L(fr)I 

-(-^)%«(i?^,)L(f.'^)L 



(59a) 
(59b) 
(59c) 
(59d) 



Here, the (3 x 6) matrices T>l and 'Dr, which enter via the diagonahzation of the (6 x 6) 
down squark mass matrix, are the exact analogues of the matrices tli and in (55). The 
partial width for the q^^ Qa^i can then be written as. 



r(r - QaZ. 



—2m„ rriy 

Ha Zi 



ML -ml - ml 

Qcr Z, 



OC; 



/ aa\ / era 



aa \ '^iJ aa^ 



xAV^ (Mj„<,m|) 



(60) 



with 



A (a;, y, z) = + + — 2xy — 2xz — 2yz. 



B. Flavour- Violating Squark Decay 

Formula (60) is very general and applies to two body decays of squarks to neutralinos for 
an arbitrary flavour structure of SSB parameters. We would expect that except, possibly, in 
the cases where we have non- vanishing values of Tq.u,d and/or Zu^d in (48), rates for flavour- 
conserving decays to neutralinos will overwhelm the corresponding flavour-violating decays. 
In many SUSY models, flavour-violating decays of squarks are thus phenomenologically 
relevant only when all flavour-conserving two-body decays are kinematically inaccessible. 
Because of the small values of down-type squark masses, this is unlikely to be the case for 
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down-type squarks. In other words, these flavour-violating decays are likely to be most 
relevant for uf^, the lightest of the up- type squarks. In many models — certainly in models 
where the SSB parameters are given by the ansatz (48) with modest values of R,, S,, W, and 
X,, and c, = 1 — the lightest charge | squark is likely to be ti, the lighter of the two squarks 
with the greatest top-squark content, and only small admixtures of ul.r and c^^ij. In models 
where the superpotential Yukawa interactions are the sole source of flavour-violation, the 
decay ti cZi has a larger rate than ti uZi because of the structure of the Kobayashi- 
Maskawa matrix. It is for this reason that this decay has received considerable attention in 
the literature [41-44]. Although other squark mass patterns are certainly possible, we will 
focus our attention on the ii cZi decay for the remainder of this section, assuming that 
the hghtest up-type squark -uf^ = ii and that the decays ti — >• tZi as well as ii — s> hWi are 
kinematically forbidden. 

Before turning to numerical details, we qualitatively estimate the various contributions 
to the effective flavour-changing coupling that causes the decay ti cZ\, beginning with 
inter-generational squark mixing. For obvious reasons, we will work in the basis that up- 
type quark Yukawas are diagonal at the weak scale. We see from the RGEs in Appendix B 
that (above all thresholds) the mixing among singlet up-type squarks occurs only via the 
Yukawa coupling matrix f„, and so vanishes (except for two-loop RGE effects) if the first 
and second generation up Yukawa couplings are neglected. Up-type singlet squarks can 
still, however, mix with doublet up-type squarks via the a-parameters. In contrast, inter- 
generation up squark mixing between doublet up squarks can occur via down-type Yukawa 
coupling matrices (which are not diagonal) even if the first two generations have vanishing 
Yukawa couplings. Thus, in models where, in the standard current basis, up squark mass 
matrices are (essentially) diagonal at the high scale, the dominant contribution to the mixing 
occurs via mixing between or t/j with c^. Within mSUGRA (or similar models), if we 
assume ~ Aof^, with ~ mo, we can estimate that the weak scale value of 

("^1)23! ~ Y|^rn2/2|K23||K33|log(MGUT/M^eak) , 

which leads to a "mixing angle" ^^^^^ ~ [vc^q)^^ /{fewx ttIq) ~ few x 10"'' [41]. With 
similar assumptions, the Ir — cl mixing angle has a comparable magnitude which, however, 
scales with the A^-parameter. Indeed, the mixing parameter that determines the decay rate 
is given to an excellent approximation in models with small flavour mixings among squarks, 
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by, 



(61) 



where the elements (^^l j (^^) (^^^"^^ J^^^ '^'^^ ~ ['^l absence 

of flavour-mixing) can be reliably estimated using the (2 x 2) submatrix of the full up-squark 
mass matrix for the mixing between and in. More generally, |e| = | (W|^)-^2 I- 

How does this "mixing contribution" compare to the "direct contribution" from the in- 
duced g-type couplings? A similar analysis to that of the previous paragraph shows that 
the dominant contribution once again arises via the and f^-type couplings. The main 
difference is that because off-diagonal elements of g arise only below the scale of the heavi- 
est sparticle, the large logarithm log(MGUT/^wGak) is absent. Moreover, the combinatorial 
factor of 8 in the previous paragraph is about unity so that these direct contributions are 
typically two orders of magnitude smaller in models where all sparticles are at the TeV scale. 
We can also check this from the magnitudes of the Yukawa couplings in Fig. 2, since the 
difference between these and the fj couplings is small. We have also checked that the value 
of the width calculated, ignoring the difference between the tilde couplings and the true 
couplings, differs from the full calculation by a few percent, as the reader may well expect. 



C. Single-step RGE integration and the stop decay rate 

The decay rate for the ti cZ\ was first estimated by Hikasa and Kobayashi [41] who 
were considering squarks of around 30 GeV that might have been accessible at the TRISTAN 
collider. They used what is essentially the equivalent of (60), but estimated the size of the off- 
diagonal elements of the up-squark mass matrix that enter (61) by a single step integration 
of the RGEs. Working within the mSUGRA framework, and using the approximation that 
Z\ ~ 7, they showed that if tree-level two body decays of t\ were kinematically forbidden, 
t\ — > cZ\ would be the dominant decay mode of ti. While their approximations^^ and 
analysis are certainly valid for m^-^ values that they considered twenty years ago, this is not 
the case for top squark masses in the range of interest today. 

Just to be sure there is no confusion, we reiterate that in models with non-vanishing values for T and Z 
matrices, this need not be the case. 

Effects of general mixing in the neutralino sector arc simple to include [42]. 
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Method 



Width (GeV) 



T{ii hWZt) 



8.6 X 10~^ 



'single-step" approximation 



41 X 



full calculation 



3.3 X 10^8 



TABLE II: Partial widths for the two- and three-body decays of the ti within the mSUGRA 
framework, with ttiq = 250 GeV, mi/2 = -250 GeV, = -930 GeV tan /? = 20 and ^ < 0, for 
which the decays ti — > hWi and ti tZ\ are kinematically forbidden. The result for the two-body 
decay is calculated using two methods, the single-step integration of the RGE (see text) commonly 
used in the literature, and the full integration of the RGEs. 

For top squarks in the TRISTAN range analysed in Ref. [41], the competing tree- level 
decays are four-body decays (we assume here that the sneutrinos are heavier than ti) of 
ti, which are both higher order in the couplings and suppressed by four-body phase space. 
However, for m^^ in the range of interest today, the three body decay tx — * bWZi may well 
be kinematically accessible and could compete with the flavour-changing two body decay. It 
is with this in mind that we are led to re- visit the rate for ti — ^ cZi, but this time integrating 
the RGEs numerically to obtain the off-diagonal elements of the up-squark mass matrix, as 
opposed to obtaining this via a single-step integration as is common practice in the literature 



To illustrate the need for integrating the RGEs we compare the rates for the flavour- 
conserving three-body decay of ti with the flavour-violating ti cZ\ decay for the mSUGRA 
point mo = 250 GeV, mi/2 = -250 GeV, = -930 GeV tan /5 = 20 and /i < (where 
we have chosen a large value of \Aq\ to obtain a light ti) in Table II. For this point, 
rrii^ ^ 181 GeV, m^^ = 102 GeV, mj^^ = 270 GeV and Wi = 197 GeV so that tree 
level two-body decays of the ti as well as three-body decays, ti — biui, to sneutrinos are 
kinematically forbidden, but both the two-body loop decay that we are considering and 
the three-body decay to bWZ^ [44] are allowed. We show the two-body decay rate, both 
for the single-step estimate as well as with the full integration of the RGEs. We see from 
Table II that the single-step approximation over-estimates the decay rate by about a factor 
of 13.6,^^ and would lead us to conclude that the branching ratio B(ti cZi) ~ 0.83, 

■^^ Since ti is mainly tn, we evaluate the decay rate using parameters at a scale equal to the lightest right- 



[43, 44]. 
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whereas the full calculation shows that B{ti hWZi) ^ 0.72, completely changing the 
qualitative picture of top squark decays! Admittedly, this striking change is because we 
are close to the kinematic boundary for the ti — > bWZi decay. We have checked, however, 
that the single-step approximation over-estimates T(ti cZi) by a factor of ~ 10 — 25 in 
mSUGRA models where the decay ti bWi is kinematically forbidden, and where the LSP 
is the lightst neutralino. This may make the four-body decay modes of ti more competitive 
than previously thought. 

D. Model dependence of r(ti cZi) 

The rate for flavour-violating squark decays will sensitively depend on whether the SSB 
parameters include genuinely new sources of flavour violation. To illustrate this, we return 
to our sample mSUGRA point from the previous section, namely mo = 200 GeV, mi/2 = 
—400 GeV, Aq = —200 GeV tan/3 = 10 and /i > 0, and include non-zero values for Tq u^d to 
obtain a variety of scenarios. We mention that the scenarios we examine are unrealistic from 
the perspective of observing flavour- violating ti decays: indeed for this sample mSUGRA 
scenario, tree-level decays of ti are accessible. Our purpose here is to understand how 
r(ti cZi) is altered (and also the quantities that it is sensitive to) as we alter the scenarios, 
to systematically allow increased non- universality and/or flavour violation in the SSB sector. 
Our results are shown in Table III, beginning with Scenario (1) which is the reference 
mSUGRA case. In mSUGRA, any dependence of the width on the matrices ^l,r{u, d) 
enters only via the KM matrix. 

1. In Scenario (2), we allow (a) non-universality by allowing T^/^) ^ but diagonal 
in our standard current basis where up- type Yukawas are diagonal at m^, while in 
(6) we allow Tq ^ and diagonal. The up-type squark matrices are (approximately) 
aligned with the up-quark Yukawa couplings so we do not expect large flavour- violating 
effects in up-squark decays in this case. Nevertheless, we do see changes of 0{1) from 
mSUGRA predictions. We emphasize that the situation would be quite different for 
flavour-violating decays of down-type squarks, since the down-squark matrix is now 

handed squark threshold. We have checked that the partial width changes by ~ 6% if instead we had 
calculated it using parameters evaluated at the highest right-squark threshold. 
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Scenario 



Width 



(1) 


mSUGRA — no dependence on specific 


^L,R{u,d) 


2.2 X 10" 


-9 


GeV 


(2a) 


Vi?(n) = YR{d) = Yl{u) = li 




T^u,D + 


3.9 X 10" 


-9 


GeV 


(26) 








1.6 X 10" 


-9 


GeV 


(3a) 


VL(n)/lL, Vn{u)=YR{d) = 


1 


Tt/,D / 


3.9 X 10" 


-9 


GeV 


(36) 






Tq/0 


2.7 X 10" 




GeV 


(4a) 


VR{d) + 11, Nr{u) = Yl{u) = 


1 


T{/,D / 


3.6 X 10" 


-9 


GeV 


(46) 








1.6 X 10" 


-9 


GeV 


(5a) 


Vr{u) / IL, VR{d) = = 


1 


Tf/,D / 


5.8 X 10" 


^3 


GeV 


(56) 








1.6 X 10" 


-9 


GeV 


(6a) 


VRiu) / IL, VR{d) / li, VUu] 




T(/,z) / 


5.8 X 10" 


-3 


GeV 


(66) 






Tq/0 


2.7 X 10" 


-5 


GeV 



TABLE III: A comparison of the two-body loop decay widths for six scenarios. For each scenario we 
list the basis used for our GUT scale inputs, with rotation matrices as specified by (51) in the text. 
In case (a) of each scenario, we take Tg = and Tj/^/) = dia^j 10000, 40000, 90000} GeV^. In case 
(6), Tu^D = and Tq = dia^l 10000, 40000, 90000} GeV^ When Tq^u,d + 0, the corresponding 
f/^jQ = 0. Here, \^L^R{u,d) are the matrices needed to transform from the current basis in 
which the matrices Tjj, T/) or Tq are diagonal at Q = Mqut to the basis where the up or the 
down quark Yukawa coupling matrices are diagonal at the weak scale. We have checked that m^^ 
is constant to within about 5% across the Table, so its contribution to the variation of the partial 
width is small. 

not aligned witli tlie corresponding Yukawa coupling matrix. 

2. This is exemplified in Scenario (36) where the mass matrix for left type up squarks 
is completely unaligned with the up Yukawa couplings. Not surprisingly, this leads 
to a very large increase in the rate for the flavour- violating decay of ti. In contrast, 
the width for this decay in Scenario (3a) coincides with that in (2a) because mfj is 
unchanged when we transform to our standard basis, so that scenarios (2a) and (3a) 
are really identical. 
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3. In Scenario (4), the rotation matrices used mean that the boundary condition for 
Scenario (4a), the Tu,d case, differs from scenario (2a) only in the boundary 
condition for m|, which, in turn, affects the running of the nig ^ matrices to a small 
extent through the RGEs. As a result, we see that the width is different from the 
other scenarios by a few percent. On the other hand, the boundary condition on iiIq 
is not dependent on V/j((i) so that scenarios {2b) and (46) are identical, and the width 
is, therefore, the same in the two cases. 

4. In Scenario (5a), the matrix mfj is unaligned with f„ resulting in a large flavour mixing 
among singlet up squarks in our standard basis, and a concomitantly large rate for 
the flavour-violating decay. It is only in this scenario that tf> — cr mixing is the 
dominant source of the flavour-violation, since in all the other scenarios that we have 
considered up to now, this mixing was suppressed by the small size of the charm quark 
Yukawa coupling. The corresponding partial width is larger than in the other scenarios 
because ti is still dominantly as can be seen from Fig. 10. (See also the discussion 
in the vicinity of this figure.) Although one might expect significant contributions 
to (59) from off-diagonal entries in the higgsino and gaugino 'Yukawa' matrices, we 



of magnitude larger than any other entry in (59). In contrast. Scenario (56) shows 
no change from (26), because the boundary condition is exactly the same in the two 
cases. 

5. Finally, the boundary condition for m|, is the only difference between Scenarios (5a) 
and (6a). Since m|, only affects the decay rate via its effect on m^g (see item 3. 
above), the change that it causes is too small to be seen in the Table for this case where 
the flavour-changing partial width is so large. Scenarios (66) and (36) coincide because 
the GUT scale boundary conditions coincide since and m|, are unit matrices, and 
so unaffected by any rotations. 

In our pedagogical examples in Table III, we considered the case with very large GUT 
scale splitting in the squark mass matrices which are likely excluded, especially when the 
squark mass matrices and the corresponding Yukawa couplings are unaligned. However, 
flavour-violation effects can be large even for very small GUT scale splitting in the squark 
mass matrices. We illustrate this for the compressed SUSY scenario, proposed by Martin [45] 




22 



term is still two orders 
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Compressed SUSY Scenario 

1.5Mj=M,=3M3, Aj=Mj, 111^=500 GeV, tanp=10, m=172 
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FIG. 11: Comparison of the partial width of the three-body tree level decay of the stop with the 
two-body loop decay in the compressed SUSY scenario discussed in the text. The dashed (black) 
line shows the partial width for the three-body decay ti — > hWZi while the other lines represent 
calculations of rates for ii cZ\ decay. For the red line in the range 10~^ we use the single-step 
RGE integration, the green line around 10~^^ is the full calculation and the upper line in the 
range 10~^ shows the enhancement as a result of taking T^^£) = (iia5{498^, 500^, 502^} GeV^ in 
the general basis using (51). 



where efficient neutralino anniliilation to top pairs via the exchange of a light squark leads 
to the observed cold dark matter relic density. We use mSUGRA-like GUT scale inputs, 
where mo = 500 GeV, = Mi, tan/5 = 10 and /i > 0, but the gaugino masses are split 
so that 1.5Mi = M2 = 3M3, with mfj^^ = diag{m'^, 500'^, 502"^} GeV^ at Q = Mqut- We 
show the partial widths for the three body and the flavour-violating two body decays of ti as 
a function of \Mi\ in Fig. 11. Over the whole of this region the two-body flavour-conserving 
decays of ti as well as its decays to sneutrinos are kinematically forbidden, but the three- 
body decay ti bWZ^ is allowed until its kinematic boundary at \Mi\ ~ 640 GeV. The 
(black) dashed line in the figure shows the partial width for this three-body decay. The 
partial width for the two-body decay, for the case of universal GUT scale masses, calculated 
using the single-step approximation is in the range 10~^ GeV and competes with the three- 
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body decay near the extreme edge of phase space. However, the corresponding result of our 
complete calculation is essentially always smaller than the width for the three-body decay, 
although it may compete with the four-body decay rate of ti at large \Mi\ [43]. 

The highest horizontal line in Fig. 11 shows the result of the complete calculation of 
r(ti cZi) for a non-universal model with Tu,d = diag {498^, 500^, 502'^} (in the general 
basis of (51)) and cu,d = 0. We see that despite the very small splitting between diagonal 
entries of the squark mass matrices at Q = Mqut, the rate for the flavour- violating two 
body decay is enhanced by three orders of magnitude, and can be competitive with the rate 
for the decay ti bWZi. We also mention that if we take this same splitting, but use 
\fi{u, d) = V/,(m) = 31, we recover the result with no splitting to within around 4%. 

Our discussion of flavour-violating squark decays in this section has revolved around the 
introduction of flavour-violation via a completely arbitrary choice of the Tujj q matrices 
in (48). It has been suggested that flavour- violation can be introduced into the theory in 
a controlled way through the minimal flavour violation (MFV) ansatz [36]. The general 
idea stems from the observation that, except for Yukawa couplings, the Lagrangian of the 
Standard Model is invariant under independent rotations in flavour space of the electroweak 
singlet up-type, singlet down-type, and the doublet quarks (along with independent rotations 
among the singlet, and also doublet, lepton flavours). If we regard the Yukawa couplings 
as spurion fields, and assume that these are the only source of flavour-violation even in 
extensions of the SM, we obtain the MFV ansatz. 

In the case of the MSSM, since m^^, and (fj^f^^)" transform the same way under the 
flavour rotations just discussed, as do iiiq and {iud^ud}^^ it is easy to see that m^^ must 
take the same form (48b) with Tu,d = 0. In contrast, for the doublet squark SSB mass 
matrix we would have, 

= ml, [11 + tjX + Urjl + ■■■] , (62) 

where the ellipses denote quartic and higher terms of the form, (f^ ^f^rf)^, (fjlfj) (f^fj) ci-nd 
(f^fj) (f*fj), etc. A similar analysis shows that. 



Auo-ti. + auiiX + /3uJlid + ---] , (63) 



where the ellipses denote higher powers of f^j^fu,^. A similar formula applies for a^. MFV 
also enters the slepton sector in an analogous manner. 
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FIG. 12: a) Variation of T{ii cZi)/m^_^ with MFV parameters Rfj (dot-dashed), (solid) 
and (dashed) that determine the form of the squark mass matrices as described in the text. 
We introduce flavour-violation into the mSUGRA model with our canonical choice of parameters 
shown on the figure. For each curve, only the one MFV parameter labelling that curve has a 
non-zero value, b) The variation of |(ti|tL)|; the ti content in the lightest top squark, which would 
equal to | cos 6t\ in the absence of any inter-generational mixing. 



It is interesting to examine how T{ii cZi) is affected if instead of choosing ad hoc values 
of Tu^D^Q as we did above, we choose these according the the MFV ansatz. Toward this end, 
we once again start with our canonical mSUGRA point, and study how this width varies 
as we switch on non-vanishing values for the MFV parameters. Since our purpose is only 
to illustrate this variation, we restrict ourselves to varying only the parameters that enter 
the squark mass matrices, and leave a,-matrices at their mSUGRA values. The variation of 
r(ti cZi)/mi^ with and is shown in Fig. 12a. The reader may be surprised by 

the large variation of the ratio (which removes the trivial growth of the width with m^-^ ) from 
the variation in Ru since we have emphasized that the SSB parameters do not introduce any 
flavour violation in this case. To understand this, we have shown the top-squark "mixing 
angle" |(ti|tL)| in frame b). Since inter-generation mixing is very small, ti is dominantly ti 
and tR, with tiny admixtures of other squarks. We see that the partial width roughly tracks 
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this mixing angle which is a reflection of the fact that for these parameters, the ti — cl 
mixing term is considerably larger than the Cl — mixing term (see the discussion near 
(61)). The fall-off in the width toward the right end of the dot-dashed curve is because the 
splitting between m?^ and the (2,2) element of M^il increases, leading to a suppression of 
e. 

Flavour violation truly enters via the SSB parameters for non-zero values of and td- 
For the tu 7^ case, the flavour-violating effects in the up- squark/ quark sector will be small 
because the up quark Yukawa coupling matrix f„ is not-very-off diagonal (since it's off- 
diagonal entries arise only from the evolution to the GUT scale) and so is approximately 
aligned with the niQ matrix at Q = Mgut-^^ The variation of the branching fraction in 
the solid curve with tu largely tracks the intra-generation t-squark mixing, just as in the Ru 
case discussed in the last paragraph. The solid curve in the left frame turns over because 
the GUT scale value of (111^)33, and hence m?^, reduces so as to suppress e. 

We have argued that the variation of the partial width for the t cZ\ decay for the 
dot-dashed and solid curves in the figure mostly tracks intra-generation t-squark mixing 
and does not derive from large flavour violation of the SSB parameters. This is sharply 
different from the result for the 7^ case, shown by the dashed curve in the figure, where 
nig receives relatively large off-diagonal contributions from the term involving the down 
Yukawa coupling matrix (which is not diagonal even at Q = rrit) in (62). Unlike the other 
two cases where the large top Yukawa coupling significantly affects the (3,3) entry of either 
M.RR or Ai-LL, in this case the ti mass and the intra-generational mixing are both essentially 
constant for the entire range of t^ in the figure, and the variation of the width is truly the 
effect of the additional flavour-violation. This is why the width continues to grow for very 
large values of \td\. The reason for the drop in the width for intermediate positive values of 
td is an accidental cancellation between the two contributions to e in (61). 

Before closing this section, we remark that the width for flavour-violating decay of ti is 
very sensitive to the parameters Ru, tu and t^. Even for values of these parameters ~ 1 — 10, 
the change in the width from its mSUGRA value is (9(10-100)%, and much larger if these 
parameters take on large values (though this may be constrained by low energy data). This 

We would, however, expect larger flavour-violation in the down type sector, since the down-type Yukawa 
matrix is not similarly aligned, and because the top Yukawa coupling is larger. 
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is a reflection of tlie sensitivity of flavour pliysics predictions to small changes in the model 
which would have little impact on the usually studied collider signals for supersymmetry. 
The flip side of this is that we must view any restriction of parameter regions from low 
energy constraints from flavour physics in proper perspective, since these will almost surely 
be sensitive to the underlying flavour structure that has a negligible effect on direct searches 
for SUSY, either at colliders or via dark matter detection experiments. 

VII. SUMMARY 

Renormalization group methods allow us to extract predictions from theories with simple 
physical principles operating at energy scales many orders of magnitude larger than the 
highest energies accessible in experiments. Due to effects arising from renormalization, these 
same simple principles lead to a complex pattern of predictions at experimentally accessible 
energies. RGEs have played a central role in the analysis of many supersymmetric models, 
generally assumed to reduce to the MSSM (quite likely augmented by right-handed neutrino 
superfields, and perhaps also additional Higgs singlets) at energy scales in between the weak 
scale and the GUT or Planck scales, where almost certainly additional new physics would 
be anticipated. 

This is the second and last of a series of two papers where we have reexamined the RGEs 
for the MSSM including threshold corrections to the one-loop RGEs (which are comparable 
to or larger than the usually included two-loop effects) along with flavour effects that arise 
below the scale of supersymmetry breaking. In Paper I [14], where we studied the scale- 
dependence of the dimensionless couplings of the MSSM, we showed that we have to extend 
the system of RGEs to include, in addition to the usually studied RGEs for gauge couplings 
and Yukawa coupling matrices, RGEs for the couplings of gauginos and higgsinos to matter 
fermions and their superpartners. Not only do these couplings evolve independently of their 
supersymmetric analogues below the scale of the highest SUSY threshold where SUSY- 
breaking effects come into play, the gaugino-fermion-sfermion coupling is also no longer 
flavour-independent (since it is no longer protected from feeling flavour-breaking effects 
below the scale of SUSY breaking) and so develops into a matrix in flavour-space. In this 
paper, we complete this program by extending the analysis of Paper 1 to the RGEs for the 
dimensionful parameters of the MSSM. 
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Toward this end, we have first adapted the RGEs for the dimensionful parameters of a 
general {i.e. non- super symmetric) gauge field theory [22] that includes interactions of two- 
component spinor and real scalar fields with one another and with gauge fields, and written 
these in a form suitable for the derivation of the RGEs for the parameters of the theory with 
four-component Dirac and Majorana spinors, together with real or complex scalar fields that 
we find more convenient for our analysis. We use non-supersymmetric methods because we 
want to include threshold corrections which, of course, break supersymmetry. The details 
of our method (including some minor corrections to the RGEs in the literature) along with 
our results for the RGEs for the dimensionful parameters for a non-supersymmetric (albeit 
not completely general) field theory are found in Sec II. Our procedure for decoupling heavy 
particles, which is essentially the same as in Paper I, is described in Sec. Ill with particular 
attention to some complications that arise when decoupling Higgs bosons, and especially 
squarks. In Sec. IV, we use these general equations to derive the RGEs for the dimensionful 
parameters of the MSSM. The complete set of RGEs is listed in Appendix B. 

In Sec. V we discuss examples of numerical solutions to the MSSM RGEs for fi and the 
dimensionful SSB parameters, focussing on flavour- violation and threshold correction effects, 
including a discussion of new technical complications that arise in the treatment of radiative 
EWSB. We have also presented the most general parametrization of high scale SSB param- 
eters valid assuming that the physics of SUSY breaking is flavour-blind. Unlike in Paper I, 
where SUSY model-dependence of the dimensionless couplings arises only via the location 
of the sparticle thresholds, the dimensionful parameters show considerable dependence on 
the underlying model. For instance, even if the gaugino mass parameters all originate in a 
common parameter mi/2, as is the case in all SUSY GUT models where the SUSY-breaking 
VEV does not also break the GUT symmetry, the resulting gaugino unification condition, 

M2 _ a2 
Ml ~ ' 

receives threshold corrections ~ 10% if the scale of SUSY scalars is around 10^ GeV, with 
gaugino and higgsino masses at the TeV scale. In Fig. 5 - Fig. 9, we show illustrative exam- 
ples of the scale dependence of trilinear scalar coupling and SSB squark mass matrices, both 
for mSUGRA as well as for non-universal models where high scale SSB parameters do not 
include a new source of flavour-violation, while in Fig. 10 we show the elements of the sin- 
glet squark matrices mfj and ml, for a scenario where a large — indeed phenomenologically 
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unrealistic — flavour-violation is introduced via GUT scale squark mass matrices. 

In Sec. VI we apply the results that we have obtained in this paper to the examination 
of flavour-violating decays of squarks. Since flavour-conserving couplings of squarks to neu- 
tralinos are presumably much larger than the corresponding flavour-violating couplings, the 
branching fraction for flavour-violating squark decays is likely to be small, unless the two 
body decays of squarks to neutralinos (or charginos) and quarks of the same generation 
are all kinematically forbidden. With this in mind, we re-visit the decay ti — > cZi that 
has received the most attention in the literature. In this connection, we flnd that within 
the mSUGRA model (where this decay has been most extensively studied) the commonly- 
used "single-step" approximation to obtain the flavour- violating iiZiC coupling typically 
over-estimates the decay rate by a factor of 10-25, and could lead to a qualitatively wrong 
picture for event-topologies from top-squark pair production for m^^ ~ 100 — 300 GeV. We 
have also examined the rate for this decay in a number of non-universal scenarios for SSB 
parameters, with and without flavour-violation. We then saw that the decay rate is sensi- 
tive to the individual matrices ^^(n, d) that enter via the diagonalization of the Yukawa 
coupling matrices, and not just to the KM matrix. Indeed, this dependence of physics on 
the separate matrices is the generic situation, while the dependence of physics on just the 
KM combination of these that we have become used to from studies within the SM or the 
mSUGRA frameworks, is true only in very special situations: see Table III.^^ We have also 
examined this decay rate in models with non-universal SSB parameters but no new source 
of flavour-violation, or where flavour violation is introduced in a controlled way via the 
so-called "minimal flavour violation" ansatz. Interestingly, as shown in Fig. 12, the decay 
rate changes by a qualitatively similar magnitude as we vary the model parameters, in both 
classes of models. The reason for this is that the variation of the width due to changes in the 
(flavour-conserving) — mixing is comparable to the true flavour- violating contributions 
to this decay rate in MFV scenarios. A determination of the left-right mixing in the t-squark 
sector (which will be difficult at the LHC [46] but possible at an e^e~ collider with sufficient 
centre-of-mass energy [47]) should allow us to readily distinguish between the scenarios. 

To sum up, in this series of two papers, we have presented the RGEs for the, in gen- 

^•^ Although this is well-known to many authors, we stress this here because there has been occasional 
confusion about this issue. For a different example, sec p. 215 of Rcf . [4]. 
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eral, complex parameters of the MSSM including one-loop threshold effects (necessary for 
two-loop accuracy) as well as flavour-mixing effects. The complete set of RGEs is listed 
in the Appendices of these two papers, and will facilitate the examination of the flavour 
phenomenology in SUSY models with arbitrary ansatze for flavour violation via the Yukawa 
coupling matrices, as well as via the SSB sector. 
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APPENDIX A: NUMERICAL INSTABILITIES ASSOCIATED WITH MATRIX 
DIAGONALIZATION 

1. The problem 

We have emphasized that in order to properly implement particle decoupling into the 
RGEs, we have to be in the mass basis of the particles being decoupled. Our procedure 
for decoupling squarks, therefore, requires us to evaluate the unitary transformation from 
the given current basis to a new current basis that coincides with the squark mass basis 
(approximated, as discussed in the main text, to be the basis in which the SSB squark mass 
squared matrices are diagonal). Below the scale Q where at least one squark has decoupled, 
we not only have the rotations ^ l,r{u, d) (unitary matrices by construction) which connect 
the current basis with the basis in which the Yukawas are diagonal at m^, but also the squark 
rotations R, which connect the current basis to the "squark mass basis" that are obtained by 
numerically diagonalizing the SSB matrices m^. R, is of course the matrix of the orthogonal 
eigenvectors of m^. If there is a degeneracy of eigenvalues, the orthogonal eigenvectors are 
not uniquely defined. This leads to a practical problem when we numerically solve for the 
eigenvectors in the case that two eigenvalues of any SSB squark mass matrix with large 
off-diagonal components are degenerate to within ~ 1%. In this case, the corresponding 
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eigenvectors, because of (system-dependent) numerical errors are not exactly orthogonal, 
and the corresponding matrix R, is not precisely unitary. 

The deviation from unitarity is very small, a part in 10^° in our case, but is nonetheless 
orders of magnitude larger than what we can tolerate when calculating the smallest off- 
diagonal elements of m^. To understand why our calculation is sensitive to this seemingly 
tiny level of noise, let us imagine what would happen if we attempted to evolve the off- 
diagonal elements of mfj from Q = Mqut in a basis where the Yukawa coupling matrices 
all have large off-diagonal elements at the GUT scale. (This is not what we actually do, but 
we could imagine doing so since we know that we are well above all SUSY and Higgs field 
thresholds where the choice of basis should be irrelevant.) In the mSUGRA framework, the 
squark mass matrices are all given by = m^'iL at Q = Mqut in any basis. Then, from 
(B17) we see that these would develop off-diagonal components ~ /ewx p x ^ x 4 x 10^ 
for the case shown in Fig. 7, where p denotes the size of the off-diagonal element of fjf*. 
In this rough estimate we have assumed that the loop factor l/(167r^) is compensated for 
by the large logarithm. In the general current basis where f„ has comparable off-diagonal 
and diagonal elements, ~ 1, and the magnitude of the off-diagonal elements of are 
(9(10^) GeV^. Rotating to our standard current basis should yield the result in Fig. 7. In 
particular, we should obtain |m^|-,^2 ~ 10"^ GeV^ because there would be large cancellations 
arising from the unitarity of Rf/ that would suppress this matrix element. If instead the 
unitarity of R holds only to a part in 10^° because of numerical errors in obtaining the 
eigenvectors, we will find that because the cancellations are not perfect all off-diagonal 
elements of will have a magnitude that is at least /ewx 100^ x 10^^° ~ fewx 10^^ GeV^, 
much larger than the magnitude of the (1,2) element in Fig. 7. We note here that the 
noise that leads to the non-unitarity of ^ l,r{u, d) matrices at the 10~^^ level is completely 
irrelevant. 

'^^ Using the g77 FORTRAN compiler with Macintosh Intel Macbook, together with the subroutine CG in the 
EISPACK collection of subroutines, we found that rIr, deviated from 11 to about one part in 10^" 
compared to a part in 10^^ for 'V L,R{u^ d^Y L,fi{u, d) for \^l.r{u, d) of the form (50). We obtain a similar 
size deviation from identity using the subroutine ZGEEV in the LAPACK collection of subroutines 
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2. The solution 



The non-unitarity of R, is only an issue when the off-diagonal entries of the squark 
mass matrix, in the basis where the Yukawas are diagonal at m^, are small compared to 
the diagonal entries. Since we, therefore, only need to consider matrices that are already 
approximately diagonal, we can associate the eigenvectors, (ei, 62, 63), with the approximate 
eigenvalues ((m^)ii, (111^)22, (1^^)33), respectively. As an illustration, let us take a case where 
the (111^)23 entry is the off-diagonal entry with the largest magnitude, and (m^)i2 the one 
with the smallest. We know the ordering quite unambiguously because above all squark 
thresholds we do not need to rotate by the matrices R, that potentially are the origin of the 
noise. We need to ensure that the (m^)i2 entry does not suffer from any numerical noise 
due to the diagonalisation. This leads us to fix ei ■ 62 = and move any non-orthogonality 
of the eigenvectors into 62 ■ 63 so that the noise moves to (111^)23, the off-diagonal element 
with the largest magnitude. To accomplish this, we slightly modify (by parts in 10^°, the 
limit of accuracy of the diagonalization routines) only the eigenvector e2 from its value as 
given by the diagonalization routines, thereby leaving ei ■ unaffected. 

To completely clarify what we have just described, although (as we have already stated) 
we do not need to rotate to the squark mass basis until we reach the highest squark threshold, 
we plot, in Fig. 13, the result for |m^|i2 obtained by the two different methods mentioned 
above, in the basis where the up- type Yukawas are diagonal at rrit, over the whole range 
Mz < Q < Mgut- The dashed (black) line shows the result where we have no rotation 
by R;7 from Q = Mgut until the highest squark threshold, beyond which we implement 
our method for ensuring that the error from the non-orthogonality of the eigenvectors of 
mfj only shows up in the (2,3) element. Indeed we see that this curve is smooth over its 
entire range of Q. The solid (red) line shows the result of carrying out the squark rotation 
without our fix of the eigenvectors over the entire range of Q. Note that there is significant 
noise all the way down to the low Q region where only some of the squarks have decoupled. 
This noise is largest at Q = Mqut where the eigenvalues of are degenerate, and settles 
down to 10^^ GeV^, not far from our estimate above. The important thing is that the 
frozen value of this element is significantly different in the two result of this 

noise, just before squark decoupling. It is for this range of Q (where the mass matrices 
that enter flavour-changing processes involving squarks will be evaluated) that we must 
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mSUGRA: mp=200 GeV, Ta^^=-400 GeV, Ag=-200 GeV, tanP=10, ^>0 

Real variables and rotations only 




le+02 le+04 le+06 Ie+08 le+10 le+12 le+14 le+16 

Q (GeV) 



FIG. 13: The scale-dependence of the (1,2) element of the Hermitian matrix for our sample 
mSUGRA point (in the basis specified by (51)) calculated using two different procedures discussed 
in the text. The dashed (black) line shows the magnitude of the smallest element, i.e. \ (m^)^^^ I = 
I (m^)^! Ii when we have used our procedure to ensure that the corresponding eigenvectors are 
orthogonal. The solid (red) line shows the same element when we do not pay attention to the 
orthogonality of the eigenvectors of m^. The lighter dot-dashed (green) line shows the magnitude 
of the (1,3) element, and provides a scale for the size of the numerical noise discussed in the text. 
The noise in this curve is too small to be visible. All elements are zero at the GUT scale. 

reduce the numerical error as far as possible. The magnitude of the (1,3) element of mfj 
is shown for comparison by the dot-dashed (green) curve. It has no visible noise because 
the corresponding eigenvalues are sufficiently split, and the corresponding eigenvectors are 
orthogonal to a very high accuracy. 

The reader will be struck by the fact that the random downward fluctuations in the solid 
curve are roughly bounded by the dashed (black) curve which shows the correct magnitude 
of the matrix element. The reason for this is that the fluctuations whose typical magnitude 
is ~ 10"^ GeV^ need to randomly fluctuate down by four orders of magnitude to even reach 
the dashed (black) line, and even more to go below, the chance for which is very small. 
Indeed it is because we have shown results for the case where the SSB squark mass matrices 
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mSUGRA: 111^=200 GeV, mj^2=-400 GeV, Aq=-200 GeV, tanP=10, n>0 



0.0015 




le+15 le+16 

Q (GeV) 

FIG. 14: Scale dependence of the magnitude of the (2, 3) entry of mfj for the same mSUGRA 
point as Fig. 13. We focus on the running at the extreme high scale, and compare the noise in the 
magnitude of just this element, both before and after fixing the orthogonality of the eigenvectors 
as described in the text. As in Fig. 13, the solid (red) line shows the result before the orthogonality 
fix while the dashed (black) line shows the result after this fix when this error has been moved to 
the (2, 3) element which now randomly fluctuates close to Mgut where the eigenvalues of are 
roughly degenerate. 

are real that we see these fluctuations go down to even the level of the dashed (black) line. 
For the more general case the chance for both the real and the imaginary part of any matrix 
element to simultaneously fluctuate downward by this large magnitude is very small, so that 
the calculated magnitude (not shown here) is always larger than 10"^ GeV^. 

After our fix of the eigenvectors, any error from the non-unitarity of R^/ is shifted to 
the largest off-diagonal element, and the only residue of the resulting noise that remains 
is in the magnitude of this element for scales close to Mqut — where the eigenvalues are 
closest — as seen in Fig. 14. At lower scales, the eigenvalues split, and the noise level (whose 
magnitude remains the same as in Fig. 13) becomes insignificant. Again, the solid (red) line 
shows the evolution of the magnitude of the (2, 3) element before fixing the eigenvectors, 
while the dashed black line shows the same thing after the orthogonality fix. We see that 
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the numerical noise has indeed moved to the (2, 3) element which now shows fluctuations, 
but only close to Mqut where the eigenvalues of mfj are roughly degenerate. 



APPENDIX B: RENORMALIZATION GROUP EQUATIONS OF DIMENSION- 
FUL PARAMETERS 



This appendix contains the RGEs with full thresholds for the dimensionful couplings of 
the MSSM with i?-parity conservation. Note that we write these RGEs in the current basis 
in which the SSB sfermion mass matrices, but not the quark Yukawa matrices, are diagonal. 
In any other basis they must be modifled to account for the rotation from this basis to the 
sfermion mass basis so that the sfermions can be properly decoupled. In this case the squark 
0q^.^s become matrices 0^^. as discussed in Sec IIIB, where further details may be found. 
The RGEs for the superpotential parameter fi and the gaugino SSB mass parameters are. 
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The following RGEs are only valid above Q = mn, where Oh = Oh = 1. We separate the 
two regimes of different Higgs boson content to simplify the resulting formulae, and to make 
explicit the parameters which remain in the theory below the heavy Higgs decoupling scale, 
Q = mu- 
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Below the scale Q = ttih, as discussed in Sec. IV B, the trilinear couplings to the doublet 
h, and the mass parameter remain in the theory, with RGEs given by, 
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With 6h = 0, RGEs for the remaining SSB scalar mass parameters take the form, 
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